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ABSTRACT 


An empirical modal of the global atomic oxygen and helium o 

distributions in the thermosphere is developed in a magnetic coordinate 

\\ 

system and compared to similar models which are expanded in geographic 
coordinates* The advantage of using magnetic coordinates is that fewer 
terms are needed to make predictions which are nearly identical to those 
which would be obtained from a geographic model with longitude and 
universal time corrections* Magnetic coordinates are more directly 


related to the major energy inputs in the polar regions than geographic 
coordinates and are more convenient to use in studies of high latitude 
energy deposition processes* This is important for comparison with 
theoretical models where the number of coordinates is limited* 


The effect of magnetic activity on the atomic oxygen distribution in 
the morning sector of the high latitude thermosphere in the auroral zone 
is also considered* A magnetic activity indicator (ML) based on an 
auroral electrojet index (AL) and the 3-hour ap index are used to relate 
the atomic oxygen density variations to magnetic activity in this region* 
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CHAPTER I 


INTRODUCTION 

V * 

( il 

l» 1 General St a f men t of tha Problem 

m ' ' ' ” ’ ’ ‘ G 

o 

Global composition and temperature measurements from mess 
spectrometer end airglow Instruments on Ogo 6 (Carlgnan and Plnkus , 

>7 ^ 

1968; Blamont and Luton, 1972) and mass spectrometers on Aeros A 
(Krankowsky et al • , 1974) Esro 4 (Trinks and von Zahn, 1975) and the 
Atmosphere Enplocnr (AE) satellites (Spencer et al., 1973; Niar et al., 
1973; Pels et al., 1973) have provided a great deal of Information, about 

the density and temperature distributions : In the thermosphere over a Wide 

0 

range of solar and geophysical conditions* These measurement's have been 
used to construct empirical models (Hedin et al*, 1974; Hedin et al., 

1977, a, b; von Zahn et al., 1977; Jacchla, 1977; Thuillier et al., 1977) 
which represent the temperature and density fields In terms of geographic 
position, altitude, and solar and magnetic activity. At high latitudes, 
Joule heating, particle inputs, and electrodynaroic forces are applied in 

very restricted areas causing large density and temperature perturbations 

c> 

over 9mall regions of the thermosphere (Hedin and Reber, 1972; Reber and 
Hedin, 1974; Taeusch and Hinton, 1975; Nlsbet and Glenar, 1977; Hisbet et 
al., 1978). In a geographic coordinate system, these energy sources 
produce large longitude effects. Towards lower latitudes, ion drag varies 

, 'O 

with the declination of the magnetic field and, therefore, induces 
longitudinal variations in the equatorial thermosphere as well (Reber et 


alt, 1973) t A combination of these polar and equatorial influences forces 
a distinction to b« made between local time and longitude# This 
diatlnction waa realised by Hedln at alt (1979) in the form of longitude 

O 

correction terms to the MSt$ (Maas Spectrometer and Incoherent Scatter) 
model and by Laux and von Zahn (1979) by aimilar corrections to the Karo 4 
model# 

An alternate approach to thia problem might consist of analysing the 
aatellite data in magnetic latitude and magnetic local time instead of 
geographic coordinates# Although the major energy input to the 
thermosphere, solar CUV heating, la organized in a geographic coordinate 
system# latitudinal and longitudinal gradients in thia energy input, 
unlike the magnetic inputs, are nail* It would thua aepo better to 
utilize a coordinate system to which the more intense localized energy 
inputs are related# If a model could be developed in magnetic coordih- f 

which did not require longitudinal correction* and Which was of comparable 

w 0 

accuracy to a geographic coordinate model developed from the same data 

bate but including longitudinal corrections, it would be very useful for 
relating energy Inputs and transport to their thermospheric consequences# 
It is the purpose of this work to develop such a model and to compare it 
to geographic representations of the neutral thermoaphere# 


1 • 2 V g«Ylous Related Studies 


U2. l V An Ovtvi»w of Thermospheric Models 


o 


Theoretical models seek to explain the observed atmospheric behavior 
consistent with the lavs of physics* Usually, a subset of the transport 
equations (continuity, momentum, and energy conservation) are solved under 
various simplifying assumptions* Boundary condition? at the turbopause 
and exobase are commonly Imposed and an isothermal-atmosphere above a 
certain height is often assumed* The difficulty in solving these 
equations has limited most theoretical treatments to one or two major 
constituents, although several authors have extended their results to 
three or more species (Vest, 1973; Harris and Mayr, 1975; Strauss and 
Christopher, 1979)* While these models are becoming increasingly 
sophisticated, usually only average conditions are treated end the 
atmospherie f tii 9 a tn|c tur e is neglected* This restriction appears to be 
easing as comparisons between theoretical predictions and actual 


observations indicate the areas in which the theoretical models ere 


deficient. 

/ n ‘ yx> 

Empirical models, on the other hand, inherently offer the advantage 
of providing better agreement with the observations since it is on the 
basis of these observations that the models are constructed* This feature 
is valuable in studies where insight into the causes of atmospheric 
variations 19 not critical, but accurate numerical predictions are 

dr . • . 

essential* Caution must be exercised, however, in the application of 
these models* Extrapolation of model predictions to conditions differing 


4 


from those under which the node! was darivsd nay give nisi sad ing results* 
although the wealth of ground-based and in situ data accumulated over the 

last decade has helped to reduce this problem* A discussion of both types 

>1 

of nodels serves to illustrate what has been done* how current problems 

C>' 

ii n 

are being approached , and what remains to be explained* 

1*2*2 Empirical Models of Temperature and Composition 

Y , ... , y “ ' // 

The early pioneering efforts of Nicolet (1961)* Harris and Priestar 
(1962) and Jacchla (1964* 196$) produced the first empirical thermospheric 
models using only satellite dreg' data* Although these models were 
landmarks in aeronomy, they were somewhat oversimplified* Since that 
time* the proliferation of aeronomy satellites and their associated 
databases in the last decade has sparked the development of several new 

O 

generations of empirical models* some of which use ground-based 
measurements in addition to in situ observations to represent the density 
end temperature fields of the earth's uppar atmosphere* I 

Alrglow measurements of the 630 nm line of atomic oxygen by the Ogo 6 

a .. - 

satellite were used to determine the global exospheric temperature 
distribution by Blamont at el* (1974)* The maximum temperature at 270 km 
was shown to lie near the summer pole during solstices where there was 
little diurnal variation. This peak migrated across the equator within 
about three weeks after the equinoxes under active solar and quiet 
geomagnetic conditions* 



4 


The earliest empirical modal baaed on measurements on several t\ 

o ’ ' , ' i ‘ V) 

Individual constituents waa davalopad by Had ip at al« (1974) from Ogo 6 Q 

mass spectrometer measurements* The data tsAse for this modal conaiatad of 

4 1 

Nj, 0, and Ha danaity masuramanta between 400 km and 600 km naar tha paak 
of tha 11-year solar cycle (27 Juno* 1969 to 13 May, 1971)* All data vara 
longitudinally (universal time) garaged, but local time variations vara 
retained* Only data taken during quiet geomagnetic periods (Ap<7 and 

ap<!2 for the same day and previous 6 hours) ware included since increased 

# ?: s " v 

magnetic nctivityrfbay have paused large fluctuations in the densities* 

/ . ; " 

Two schemes vara used to model the composition; in one, the densities ware 
extrapolated to 120 km, using exospheric temperatures and a temperature 


w 


•l 


/ 


gradyiunt parameter (s-0.0215 km ) determined from the N 0 measurements and 

"••••; /S ... ** 

based on a Bates (1959) temperature profile having constant lower boundary 


/ 


(120 km) values* In the other, the densities were extrapolated to 450 km 

/> 

aasumlng an isothermal regime* A least squares fit of spherical harmonics 
to the data was made in both cases* These functions were used since it 
waa felt that they were approximate eigenfunctions in the thermosphere and 
since they were complete. That is, only a small number of terms should be 

needed, hut as much detail as desired could be obtained simply by adding 

/ ' // 

more terms. The theoretical advantages of spherical harmonics are 
discussed more completely by Mayr and Volland (1971) and Volland and Mayr 
(1972)* The model contained only time independent, solar and magnetic 
activity, annual, semiannual, diurnal, semidiurnal, and terdlurnal terms* 

’ O' 

The variation included nonlinear components, hut because the range 

was small, Ap variations were linearly approximated and a 6-hour lag was 


l 


// 






- . . 
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used In accordance with Jacehia (1971) . 

Tha accuracy of tha temperature representation was estimated to b a 
+50 K subject to tha validity of tha boundary conditions* Ilia accuracy of 
tha 0 and Ha value* at 120 km were dependant on tha assumption of 
diffualva equilibrium as vail as on tha infarrad tamparaturaa* Tha 
uncartainty in tha actual measurements of tha densities, up to 252 for 
and 0 and up to 502 for Ha, not including a 10-152 calibration 
uncertainty, also limited tha modal's accuracy* 

A diurnal maximum near 1000 LT was found for Ha while Ng peaked near 
1600 LT. Winter maxima of Ha and of tha 0 to Nj ratio vara also observed 
in addition to a 400 K summer to winter exospheric temperature difference. 
The Inferred temperatures appeared to be about 72 higher than 
corresponding incoherent scatter measurements (Salah and Evans, 1973) and 
the diurnal variations compared favorably with the exception of the winter 
morning values trtiere Ng data were sparse* 

’’ 'o 

Kasprzak and Newton (1976) compared predictions of the Ogo 6 model 
(Had In at al«, 1974) to mass spectrometer measurements from San Marco 3-> 
They found that the Nj values agreed within +402 while total oxygen and 
mass density were within +302; however, poor agreement between helium 
values was found and attributed to extrapolation of the Ogo 6 data to the 
lower altitudes which were covered by San Marco 3* The Ogo 6 model showed 
a diurnal maximum of He from 1000 to 1100 LT at 450 km while San Marco 3 
data displayed a maximum from 0700 to 0800 LT between 220 km and 280 km* 




m 


a 


men, 


7 

0 

A modal similar to Ogo 6 (Had in at al., 1974) was presented by von 
Zahn at al* (1977) on the basis of Esro 4 composition measurements* This 
modal used data over the 240-320 km altitude range and covered the lover 
activity portion of the solar cycle* Exospheric temperatures were 
Inferred from both molecular nitrogen and argon* These averaged 5-92 
lover than those predicted by the Ogo 6 model for an F 10. y of 120 and were 
supported by more recent findings (Hedln et al*, 1975; Mauersberger et 
al*, 1976; Thuilliar etal*, 1977)* The summer to winter temperature 
difference was also lover for the Esro 4 (315 K) than for the 0$o 6 model 
(415 K). The seasonal a N n^, diurnal variations of N, and Ar In the Eero 4 

j) 2 

model were similar, supporting the assumption that both reflect the 
thermal structure in the upper atmosphere. The seasonal variation of Q 
appeared be dominated by a semiannual effect and the diurnal amplitude 
was found to be small* The He distribution showed a noticeable winter 
bulge of greater than 40 tjaes the summer density and the semiannual 
component was found to be much larger than the diurnal variation* 

An extension to the Ogo 6 model was made by Hedln et al* (1977a, b) 
using mass spectrometer data from five satellites (AE-B, Ogo 6, San Marco 
3, Aeros A, and AE-C) and Incoherent scatter measurements from four ground 
stations (Arecibo, Jlcamarca, Millstone Hill and St* Santln) * The average 
exospheric temperature for this MSIS (mass spectrometer and Incoherent 

scatter) model was 1041 K for an of 150* This was 56 It lower than 

§ . 

the Ogo 6 model and 35 It higher than the model of Jacchla (1971). The 
uncertainty of this value ves estimated to be 15-20 K from comparisons 

between Incoherent scatter temperature and N, density data. The low 

0 


altitude data also agreed wall with rockat data (Offar&snn, 1974 ) but 

tha higher altitude Ng values differed markedly from the Jacchla (1971) 

o . C. 

values* The othe!} apeciee represented (lie, Ar , 0 g and H) and , in 
particular* 0* did not agrae so Well with low altitude rocket data 
(Offeinuann, 1974)* possibly due to surface reactions of 0 in the rocket 
data and incomplete geographical coverage of He* The high altitude He 
variations were nearly th/l same as found by Keating at al* (1974) on the 
basis of satellite drag analysis and the seasonal variations of atomic 

v 

oxygen agreed well with the results of Mauersberger et al* (1976)* 


Thuillier etal* (1977) extensively analyzed the 630 nm air glow data 
from Ogo 6 and presented a global empirical model of the exospheric 
temperature derived from this data. The polar heating patterns seen by 
Blamont et al* (1974) were still present after the database of the model 
had been edited and these patterns occurred during both solstices and 
equinoxes* They were particularly strong during magnetically disturbed 
periods* The alrglow derived temperatures were in fair agreement with 
Incoherent scatter results (Salah et al*, 1975) and with temperatures 
Inferred from Ng density data (Hedln et al*, 1974) and Ar and N, densities 
(Chandra and Spencer, 1975). 

A comparison of exospheric temperatures inferred from total density, 
N, density, and 630 nm alrglow data made by Nlsbet et al. (1977) showed 
that temperatures inferred from Mg densities were in much better agreement 

with alrglow temperatures than they were with those Inferred from total 

0 ... 

densities* This was especially true at high latitudes where the 


L, 


’ \ 

atmosphere is definitely not in diffusive equilibrium* However, the mess 
of is closer to the mean molecular mass than any other constituent and 

N 2 is the major constituent in the region where eddy diffusion and 
vertical diffusion ere both important and it is a batter Indicator of 
temperature than any other species, including the sum of all gases, under 
t|4 assumption of diffusive equilibrium* The H 2 density inferred 
emperatures did, however, underestimate the influence of magnetic 
activity on the temperature at low latitudes while overestimating this 
effect at high latitudes using CIRA (1972) boundary conditions at 120 km 


and the polar cap distributions of N, and the temperature were rather 


O 




different. 


The effect of the large scale circulation induced by a magnetic storm 

in February, 1974 on the neutral composition was investigated by Mayr and 

Hedin (1977) utilizing the model of Mayr and Volland (1973a), High 

latitude depletions of 0 and He and enhancements in Ar as well as low 

latitude Increases in 0 and He and decreases in Ar were predicted and 

substantially validated by AE-<3 neutral density measurements. An energy 
~3 -2 

input of 1*7 X 10 ,W m above 120 km was needed to bring the theory into 
agreement with the composition data. The annual changes in the neutral 
density distribution were said to have an Important influence on the 
dynamics of magnetic storms* It should be noted that this treatment was 
based on a model having no local time dependence either in the energy 
input or In thie neutral density. 
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Jacchis at al* (1977) used composition measurements from Earo 4 to 
ccaata a global thermospheric nodal for nag na tic. ally disturbed conditions* 
Thay showed that a variation of the Haight of the horaopeuae would 
reproduce the density changes observed during disturbances* though noting 
that rapid fluctuations of this level night not be likely* They also 
found that the differences between density values obtained during 
transient disturbances and sustained nagnetic scorns could be reproduced 
by the nonlinear character of the relationship between geomagnetic heating 
and the variation in homopause height* The sudden heating gave rise to an 
'equatorial wave' in addition to the usual thermal effects* 


Barlier et al* (1979) compared several thermospheric empirical nlobal 

\\ ,j 

nodels (MSIS^\\Hedin et al*, 1977a, b; Esro 4, von Zahn et al*, 1977; J77, 

. v 

Jacchia, 1977; &TM, Barlier et al*, 1973b) under various geophysical 
conditions* Each model performed relatively well overall, but none of the 
models explained the variations present in the data completely* In 
particular, diurnal phase and amplitude estimates were noticeably 
different below 200 km* Also, major discrepancies between He deta In 
areas of low He concentrations from satellite drag analysis and from mass 


spectrometer measurements could not be reconciled* 


Hedin et al* (1979) and Laux and von Zahn (1979) made improvements to 
the MSIS and Esro 4 models, respectively, by adding longitude corrections 
to these models which reflected the offset of the geomagnetic poles from 
the geographic poles* These corrections represented the gross effects of 
magnetic coordinates on the composition and density, but the role of « 


increasing magnetic activity was still not satisfactorily delineated* 

Tha annual variations of savaral neutral constituents (N,, 0, Us, At, 
N) vats outlined by Kbhnlsin at ai . (1979) using Aatos A and B data taken 
near 0400 IT and 1600 LT* Agreement with MSIS and Esro 4 predictions was 
quite good for He although a large scatter (+401) was present in some of 
the measurements* 'the general behavior of N agreed fairly well with 
Engebretson et al. (1977) although Aoroa maaurementa displayed a polar 
structure which was quite different* The Aeroa, Euro 4, and MSIS models 
predicted nearly identical values for 0 except where the measurements were 
extrapolated to 120 km* The Aeros date on Ar and Nj showed strong 
variations with magnetic activity near the poles as was predicted by Esro 
4 and MSIS, hut the equatorial behavior of Ar waa quite different then 
MSIS estimated* 

£ 

1*2*3 Theoretical Models 

O 

O j 0 

1* 2* 3* 1 Neutral Composition and Temperature 

The seasonal behavior of the neutral thermosphere averaged over daily 
variations was examined by Jfayr and Vo Hand (1972a) Who outlined the 
latitudinal structure of the annual and semiannual components by applying 
small perturbations to the conservation equations and observing the 
effects of different energy inputs. EUV heating Was used as art energy 

' . m - - _ - 

source for the annual component and a term simulating auroral heating, 
which maximized at the polea, was added to this in order to study the 
semiannual variations* Although photochemistry was neglected, clear 

° O O'. 
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annual and semiannual trends vara found* The annual aolar haafc Input 
Imbalance induced a circulation system with wind* bloving from the summer 
to th4 winter hemisphere in tha uppar thamoaphara* Thia first appeared 
to damp tha tamporatura variation strongly; however, in ordar to maintain 
continuity* the lightar minor constituents 0 and He began to diffuse 
upward through tha major gas molecular nitrogen in tha summer haniaphara 
and became depleted at lovar altitudes where the large molecular nitrogen 

concentration adted as a barrier to diffusion thus decreasing the pressure 

' f ' 

variation and, consequently^ the horizontal wind velocity. Thia in turn 
reduced the adiabatic heat transfer from the preferentially heated auramer 
hemisphere to the cooler winter hemisphere allowing the maintenance of a 
larger temperature imbalance* The semiannual affect appeared as reduced 0 
densities below 450 km and reduced He densltlee at all therraoepheric 
heights near both poles corresponding to tiff |-aak in heat input 
there; however, the equatorial concentrations of these constituents 
Increased in response to the semiannual circulation cells from the poles 
to the equator* Surprisingly enough* the annual mass density variation 

o . 

was found to be nearly constant with altitude and the small semiannual 
variation almost latitude Independent; these results were due to the 
circulsti< 3 n> : ,and the resulting diffusion of atomic oxygen and helium 

. C. 

a • • . 

through molecular nitrogen* This explained for the first time why 
satellite drag measurements of total density showed little change in the 
annual variation with height while mess spectrometer data indicated the 
opposite* 
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The diurnal phase dlfferanc* between the temperature and the 0 » 

O . ... *\ 

density was Investigated by Mayr and Volland (1972b) uairg a 0 
two-dimensional time-dependent model having solar heating as an energy 
sourcei the momentum and continuity equations were solved assuming no 
variation In n(0) at the lower boundary of 90 kra and diffusive equilibrium 
for atonic oxygen above 250 km. Wind-induced variation* of n(Q) were 
found to control the atonic oxygan distribution below 200 km whereas 
thermal expansion became more Important abova 300 km. Tho time of the 
diurnal maximum of atomic oxygen of 1030 LT was nearly equal to that of 
the wind at low altitudes. Increasing slowly towards that of the temper- 
ature (1600 LT) throughout the upper thermosphere. The authors noted ' 

C> 

that a major shortcoming of this model was that energy Inputs from mag- 
netic storms, which were not Included, could couple with the lower 
thermosphere and possibly reduce the density-temperature phase differ- 
ence significantly. Also, the effects of semidiurnal and terdiurnal 
harmonics, which could be important, were not considered. 

/ 

This study was extended by Mayr at al. (1973), who used a 'quasi' 
three-dimensional model and considered diurnal, semidiurnal, and 
terdiurnal harmonics in the energy input term. These harmonics manifested 
a latitudinal dependence and peaked at the equator. Nonlinear coupling 
between the diurnal ion drag term and the diurnal horizontal wind 
component was also included. The 1600 LT temperature maximum and the mass 
density peak near 1400 LT were shown to be consistent with the theory. 

The semidiurnal and, more importantly, the terdiurnal harmonic caused a 
shift in the temperature maximum at high altitudes to later local times in 
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clot's agreement with radar baqksqatter observations (Carru at al*» 

1967? Niebet# 1967; Waldteufel and McClure, 1969? Mshejen, 1969? Salah and 

' ') 

Evant, 1973)* Lower In the thermosphere, energy redistribution raaulting 
Iron heat advactlon tandad to ahi£t tha temperature paak towards noon* 
Diffusion of 0 and Ha, especially from anargy inputs occurring below 120 
ka, greatly influenced tha phase of tha mss density above 250 ka» 
shifting it towards earlier hours* 

\ The effects of upper atmospheric flow on the minor gasos helium and 
argon were outlined by Rebar and Hays (1973) by combining the momentum and 
continuity equations to derive density profiles for these gases* Modal 
valuas ware uaad for the temperature and major gas (N^) distributions and 
lower boundary densities were assumed for He and Ar* An expression 
developed by Hodgds and Johnson (i960) was used to determine the upper 
boundary flux of He* Ap arbitrary vertical wind profile which Induced / 
summer to winter transequatotial flow was chosen as the thermospheric wind 
field* This hemispherical flux caused a significant increase in the 
winter to summer helium density ratio* Exospheric transport lessened this 
value, particularly near the peak of the linear solar cycle when the 
exospheric temperature was high, but a pole to pole ratio of about 10 was 

' ■ »i 

■till maintained* Equatorial horizontal winds of 100-200 m sec * at 200 
km were found to be consistent with the He densities measured by Ogo 6* 

The effects of the wind field on Ar were quit) different. The Ar density 
increased in the summer hemisphere by almost a factor of 4 over winter 

hemisphere values and exospheric flow was found negligible since the Ar 

\ ' ” 

density decreased so rapidly with height* No measurements of Ar were 
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available to compare with this estimate at that time* 

The diurnal behavior of the thermosphere war. examined by Harris and 
Mayr (1975) by extending the previous studies of Mayr and Vo Hand (1972a, 
b; 1973a) involving nonlinear processes and 'mode coupling'* Mpde 

coupling was believed to occur between 100 km and 200 km where the 

o 

functional form of the solutions to the transport equations changes 

O ° ' 

radically* In this transition region, solutions could only be expressed 
by a combination of both low altitude and high altitude eigenfunctions* 

The effects of viscosity, Ion drag, diffusion, advection and heat 
conduction were also Included In the analysis* The results showed a 
decrease -in the diurnal variations of temperature and Nj density at all 
altitudes with respect to one-dimensional models* The 0 variation was 
more complicated; above 200 km the diurnal amplitude decreased, 0 but below 
200 km a substantial Increase (due to wind-induced diffusion) Was 
observed* The phases of the temperature, density, and 0 density 
decreased at nearly all heights; in particular, the phase of 0 below 200 
km was shifted by more than 12 hours towards morning* The unique 
amplitude and phase behavior of 0 were direct results of transport* 

Further, these results were what Harris and Priester (1962) had attempted 

ij: ' 

to duplicate by Invoking a ' secoqi heat source'* 

Strauss et al. (1977) constructed a three-dimensional dynamical model 

Jj 

of helium in the upper atmosphere which Included the effects of exospheric 
transport, molecular and eddy diffusion, and solar activity. A background 
gas consisting of 0, 0^, and N 2 was assumed and represented by the CIRA 
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(1972) model, and momentum feedback by He on this gas was considered 

<? 

negligible* The empirical representation of the background gas density 
and the temperature was considered more realistic then the solution of the 
energy equation due to uncertainties in the solar flux estimates* The 
only equations requiring solutions were, therefore, momentum and 
continuity for Ha* An upper boundary condition for helium was taken from 
Hodges and Johnson (1968)* The low altitude diurnal maximum in helium was 
estimated to lie between 0600 and 0800 LT» Although this phase Increased 
somewhat with height, it was still earlier than empirical models 
predicted* Lower spatial resolution in the theoretical model was 
suggested as an explanation* Tha pole to pole helium ratio was calculated 
for aevkfal values and was found to increase as the ^ index 

decreased* A decrease in exospheric transport seemed to be the reason for 
this as Rebar and Hays (1973) had found* Smaller effects on this 'bulge 

■ vV. 

ratio' were observed by varying K, the oddy diffusion coefficient* The 
ratio increased from 7*2 to 9*0 at 4^P km as R was increased from 3 x 10 
to 10 3 m“ s'* 1 for an F^q^ of 144 x 10" ta ^ W m"** He”*. These results were 
generally in agreement with empirical models and quantitative differences 
were attributed to insufficiences of the CXRA (1972) parameters used* The 
authors suggested that the inclusion of a globally varying eddy 
diffusivlty could have Significant effects on the predicted distribution 
of the minor gases* 

Strauss (1977) estimated the global distribution of argon by the same 
method used to determine the distribution of helium, but without 
exospheric transport* Calculations were made for both hydrostatic and 
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dynamic conditions. For dynamic conditions, the dlurnsl maximum of argon 
appaared latar In the day (1400 LT) and the diurnal 'bulge ratio* ffi$a 
larger (U.5 at 320 km). The model compared favorably in the equatorial 
regions with San llarco 3 data (Newton et al«, 1975) in determining the 

diurnal phase; however, the diurnal amplitude waa overestimated and an 

j) s 

inadequate repreaentation by the CIRA (1972) model, which was used to 
calculate the vertical diffusion velocity of argon, waa suggested as the 
reason for this. Comparison with the MS IS model (Hedin et al., 1977a, b) 
and with data from Earo 4 (von Zahn et al*, 1/973) again revealed 
overestimation of the diurnal variation. 


Strauss and Christopher (1979) investigated the Influence of 
transport on the global atomic oxygen, molecular oxygen and molecular 

nitrogen distributions in the thermosphere. Using 0, 0^, and N 0 coupled 

o 

by collisions and neglecting photochemistry, calculations for hydrostatic 
and steady state dynamic conditions were made by solving the mass and 
momentum conservation equations* Departures from diffusive equilibrium 
for Og and N fJ below 200 km Were alight, but atomic oxygen showed a 
substantial variation In this region. All constituents approached 
diffusive equilibrium towards higher altitudes. A low altitude minimum of 
the 0 concentration was found at both poles for equinox and at the summer 
pole for solstice. The local summer minimum gave way to a maximum in 0 at 

II 

higher altitudes (>375kra). These effects were consistent with 


wind-induced diffusion in the lower thermosphere succumbing to thermal 

0 . . ■ .. 

expansion with increasing height. 


s 



i 

s 



1 


% 


t! 

V ^ 


0 


18 


1.2. 3. 2. Atmospheric Motion 


Theoretical examination* of the thermospheric circulation syatem 
have traditionally been based on the solution to the equation of motion 
for the neutral gas using empirical temperature and pressure distribu- 
tions and assuming diffusive equilibrium. Although this assumption has 
been found unrealistic below 150 km (Mayr and Volland , 1972a, b; Mayr 
et al. , 1973), it has, nevertheless enabled investigators to produce 




important and useful models of the teg 
investigated by satellites. 


bn of the thermosphere normally 


One of the earlier attempts to represent the thermospheric wind field 
was made by Kohl and Kins (1967) baaed on a static atmospheric model 
(Jacchia, 1965) used to estimate pressure gradients resulting from 
non-uniform solar heating in the upper atmosphere. "These estimates of the 
'driving force' were applied to the linearized horizontal equation of 
motion for the neutral gas which included ion dtag, viscous, Coriolis and 
linear inertial terms. Vertical profiles of the horizontal wind velocity 


in the absence of electric fields were thereby obtained. The calculations 


1 

were performed for equinox for peak electron densities of 10 * ra , 


|| 

roughly corresponding to daytime and 3 x 10 m , approximately those of 


nighttime. The wind velocity vector at 300 km corresponding to high 
electron densities essentially followed the pressure gradient due to the 
strong influence of Ion drag* The Coriolis and inertial terras dominated 

the circulation pattern for the lower electron density (lower ion drag) 

o y ' 

conditions. In this case, the high latitude ^nd direction was distinctly 
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shifted towards the east while the low latitude motions were directed 
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westward* The authors cautioned that these results were subject to the 
rather restrictive assumptions of a sinusoidal driving force* fixed 
electron and ion densities* a aingla viscosity coefficient profile, 
coincidence of geographic and magnetic polee* and tha neglact of all 
forces other than that^off the pressure gradient (e.g. electric fiolde). 

A similar study was undertaken by Ce taler (1967) using a modified 
pressure distribution (Jacchia and Slowey* 1967) as the driving force* 

The resulting wind patterns did not vary substantially from those found by 
Kohl and King (1967) although the vind velocities differed* Variations In 
ion drag appeared to have much more influence on the meridional wind 

o 

Component than changes in molecular viscosity* 


Atmospheric models for both neutral and ionized partlc 


Marts' Mas 


(dttilchiag.. 1965) Nisbet, 1970) formed the basis Cor the global wind system 
study completed by Blum and Harris (1975a* b)« The Navier-Stokas 
equations including nonlinear effects were solved to obtain horizontal 
winds* although the nonllnearltles were subsequently shown to be of 
importance only in the equatorial region* In contrast to previous work, 
the effects of electric fields were considered, but the simplification of 
coincident geographic and magnetic poles was retained* The pressure 
gradient and ion drag terms appeared to have a dominant influence on the 
solutions, but significant variations from Coriolis, viscosity, and 
inertial terms at various geographical positions indicated that 
simplification of the equations of motion was only valid for fixed 
altitudes or geographic locations* Comparison of theoretical predictions 
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with wind velocities determined fr om incoherent scatter measurements vara 
generally favorable* the empirical findings of Harper (1971) at Areclbo 
(18°M) indicated smaller aqua Cot ward velocities than predicted and a 
diurnal phaaa about ona hour earlier than the thaoraticai estimate* The r 
observations of Evans (1971) at Millstone Hill (42*6°H) also ahowad a 
smaller equatorward average but tha lower ordar harmonica of tha local 
time variation were similar to tha thaoraticai curvaa. Tha meridional 
winds deduced by Aaayanc and Vaaseur (1972) for equinox also had a smaller 
equatorward maximum than the theoretical aa Ornate a* This discrepancy was 
not so largo for solstice and was thought to hava bean due to tha lower 
ion densities used in the theoretical treatment* 

A similar approach was taken by Craekmora at *1. (1975) who used the 
Jacchla (1971) modal for tha temporatura and neutral densities and the 
Ching and Chiu (1973) model Cor tha ion concentrations* Modes up to the 
4th harmonic vara included in the treatment and the results were somewhat 
similar to those of Blum and Harris (1975a, b) when the differing ion 
densities assumed were reconciled* large scale effects appeared to be 
illustrated reasonably well, but the resolution needed to elucidate the 
/Smaller scale features was lacking* The significant contribution of 
vertical motion to the energy balance at high latitudes was stressed as an 
important area requiring future work. 

The influence of a convection electric field in the polar regions on 
ion drifts and neutral winds wet investigated by Maeda (1976). A dawn to 

dusk electric field similar to that observed by Ogo 6 (Gurnett, 1972) and 
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two height Independent latitude profiles of the Ion denelty were used to 
solve the horizontal equations of motion for the neutral and ionised gases 
between 110 km and 200 km* The neutral wind field of eleatrlc field 
origin was found to depend strongly on the ion density profile and the 
wind speed was found to be comparable to the thermally driven wind speed 
(Kohl and King* 1967) of greater than 100 m sec ^ at 200 km* 


Roble et el* (1977) analysed the seasonal and solar activity induced 
variations of the zonal mean circulation for quiet geomagnetic conditions* 
Extending the model of Dickinson et al* (1975* 1977) to low solar 
activity* the authors found that the mess flow above 100 km at solstice 
from the summer to the winter hemisphere was much larger near solar 
maximum than near solar minimum* The equinoctial transition of the 

O 

average in terhemi spheric flow occurred within 2 weeks of equinox and this 

0 

delay diminished substantially with increasing geomagnetic activity* 
Although the flow was driven mainly by solar UV and EUV sources* 
modulation of these sources by high latitude effects which were themselves 

related to geomagnetic activity appeared likely due to the variation in 

// 

this transition period* Since the magnitude of the high latitude sources 
needed to be adjusted over the solar cycle to obtain agreement with 
measured temperatures and winds* the energy input In these regions was 
expected to vary with solar activity* 
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1.2.4 Small Seal* Phenomena 

1. 2.4*1 Temperature and Penalty Behavior Purina Periods 
of Low Magnetic Activity 

The special and temporal variations of the neutral thermosphere which 
are highly localized have not yet been Included In global empirical 
nodela. The large number of terms needed to adequately represent these 
perturbations In an analytic fashion do not seem justified at this time 
because of Insufficient physical understanding and limited databases. 

Only specific Inquiries Into regionalized phenomena can give an accurate 
picture of the fine structure In the thermosphere at the present time. 

Mid- to hl'gh-latltude heating sources were observed by Reber and 

jj 

Hedln (1974) by analyzing two years of neutral density measurements from 
Ogo 6. The molecular nitrogen concentration near 450 km was found to 
Increase 1 the helium density Co decrease, and the atomic oxygen 
concentration to remain fairly stable in the region of heating. These 
results appeared consistently during magnetically quiet periods, 
suggesting that the heating was a permanent feature of the neutral 
thermosphere. Corotation of this heat source with the geomagnetic pole 
Indicated that some type of geomagnetic control was being exerted although 
no correlation with the cusp region could be found. 

Taeusch and Hinton (1975) related N, density variations to 
electrodynamic and particle heating in the quiet polar thermosphere. Mass 
spectrometer measurements of the Nj density from Ogo 6 were normalized to 
450 km and then divided by the corresponding Ogo 6 model densities (Hedln 
et al . , 1974) In order to estimate- differences from the average 
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concentration as Well as variability in these differences* The behavior 
of N 0 with respect to ‘invariant latitude (ILAT) and magnetic local time 
(ML?) was then outlined* The regions having the greateat differences 
between measurementa and model estimates were the cusp area and the 
0000-0300 ML? sector* Although the position of the cusp was quite 
variable, it was usually located near 80° ILAT and 1200 MLT. In addition 
to large differences from average values* the densities in the cusp were 
found to be quite variable* The densities in the 0000-0300 MLT region did 
not display this variability* This suggested that the energy inputs to 
the cusp, quite possibly the particle precipitation identified by 
Gustafson (1972), were quite variable while the heat sources in the 
0000-0300 MLT region, which ate more likely electric fields over a region 
where particle precipitation has enhanced the conductivity, ware very 
persistent* 

Wavelike structure in Ar, and He was seen by Reber et al* (197S) 
over several individual orbits of Atmosphere Explorer C (AE-C) ti The 
relative amplitude of the Ar waves was about twice that of the 
disturbances, but with the same phase* He was usually 130° out of phase 
with Ng and had about one-half the amplitude* The data examined, 
however, were too limited to determine the nature of the disturbances/! 

I » 

The first experimental evidence of transport by the diurnal tide was 
presented by Hedln et al* (1978) using AE-E data* A large phase shift in 

fj \ . 

0 and He from high to low altitudes, predicted by Mayr and Harris (1977) 
was substantially validated by these observations although marked 
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statistical deviations in the data mads it difficult to svaluats the 

* 

agreement between theory and observations quantitatively* the region 
above 200 km showed no departures from diffusive equilibrium while these 
were marked in the low altitude regime* A description in terms of eddy 
diffusion models, which predicted composition variations only below 120 km 
and diffusive equilibrium above, was found unsatisfactory, since the 
predicted Increases of 0 and He amplitudes from 160 km to 120 km were not 
observed . 

Temperature variations in the polar thermosphere over the 11-year 
aolar cycle were studied by Chanin and Tulinov (1979). Sodium clouds 
between 120 km and 180 km from rocket flights were used to estimate 
neutral temperatures* Values obtained near 165 km for the solar minimum 
were found to be about 500 K larger than those for the solar maximum* A 
possible explanation of this was increased turbulence in the lower 
thermosphere and therefore, larger heat flow to the mesosphere during the 
aolar maximum (Roper, 1973)* This flow would cause a. higher altitude 
meaopause at the peak of the solar cycle and, therefore, lower 
temperatures near 120 km* This explanation would be consistent with the 
observed higher exospheric temperatures of the aolar maximum since the 
temperature gradients Would be larger in this case* 


l*2*4*2 Temperature and Density Behavior Purina Periods 

Analyst* of orbital drag measurements during geomagnetic atoms by 
Jacchla (1959) first auggaatad that the neutral thermospheric density 
Increases during these disturbances* The tine lag of the response of the 
composition to heating as well as the details of the response continue to 

be a topic of considerable Interest* .Although the theory is still 
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developing* a summary of the advances made is In order* 

A theoretical investigation of the density response to auroral 
heating was launched by Vo Hand and Mayt (1971)* A heat input resembling 
an impulse function in latitude and exponentially decaying in time was 
positioned in the nighttime auroral zones {+ 65°LAT> and used to drive a 
three-dimensional model of the thermospheric density based on a spherical 
harmonic expansion to see which harmonic components were effectively 
damped out and which propagated to lower latitudes* Although an 
inordinate number of harmonics would have been needed to accurately 
portray the heat Input* it was shown that the thermosphere effectively 
damped all but the four lowest components after the first hour of the 

storm* The time response of the maximum density change varied from 2* 5 

C' 

hours at the pole to 6 hours at the equator and the amplitude of this 
variation at the pole was twice that at the equator* The local time 
behavior of the density changes was found to depend on latitude, 
maximizing between about 0000 LT and 0600 LT. 
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Neutral densities from the Ogo 6 mass spectrometer were examined by 
Taeusch et al# (1971) to determine what variations of composition occurred 
during the geomagnetic disturbances in the period 27 September to 3 
October* 1969. All data were normalized to 500 km and divided by quiet 
time values before comparisons were made. The ratio of the N 2 density 

under disturbed conditions to the N 2 density under quiet conditions 

, ■ !/ 

showed the largest value in the auroral zones and the 0 ratio increased 
here as well* but not by so much. The equatorial variations were small 
during the storms observed. The N 2 and 0 density variations were used to 
infer neutral temperature changes by invoking Jacchia (1971) model 
temperatures corresponding to the neutral density values. Temperature 


i of 450 K ot ysii }!- whila iafigiiili 


of only 325 K were inferred from 0 variations. The smaller increases for 
0 were attributed to changes in the assumed lower boundary conditions. 

That is* N 0 was expected to be more nearly in diffusive equilibrium at low 
altitudes than 0. In contrast to the high latitude temperature increases* 
equatorial enhancements of only 25 K for N 0 and 45 K for 0 were 

m 

calculated. These values suggested upwelling in the polar regions and 
subsidence near the equator as predicted by Mayr and Volland (1970) and 
others* The density changes above 60° geographic latitude appeared to 

V y ;; , 

occur less than one hour after the onset of a storm as measured by the 
geomagnetic field variation* The largest changes occurred near 0600 and 
1000 MLT which suggested geomagnetic control of the heat source(s). 


A decrease in the lower boundary 0 to ratio and changes in 
thermospheric circulation during magnetic storms were predicted by Chandra 
and Stubbe (1971) upon solving the equations of notion and heat conduction 




for the Ion* and neutrals and the ion continuity equation* The height of 
the turbopause was expected to vary and any temperaturea derived ualng 
fixed boundary condltiona were* therefore* considered to be inaccurate* 

The ion densities were shown to be sensitive to the 0 to ratio and it 
was suggested that ionospheric parameters be used in modeling since they 
were easier to measure than neutral ones. 

Hayv and Volland ( 1972c) showed that the circulation system set. up by 
a magnetic storm would cause low altituda depletions in n(0) in the polar 
regions due to diffusion and high altitude enhancements in n^)* but not 
in n(Q), at high latitudes* The model of Mayr and Volland (1971) was used 
for the calculations and two auroral heat sources were assumed to excite 
thermospheric circulation* One source maximized near 100 km* simulating 
auroral particle precipitation and the other near 150 km* which was 
assumed to be similar to Joule heating by electric fields* A constant 0 
density at 90 km and no diffusion above 250 km* where 0 becomes the major 
species* were assumed and solutions to the coupled continuity and momentum 
equations were obtained using up to the second Legendre polynomial for the 
heat input* temperature and wind fields, and the diffusive flux* Only the 
Joule heat source was found capable of reproducing the densities measured 
near 400 km during a magnetic disturbance (Taeusch et al** 1971)* The 0 
density variation was separated into temperature and wind components* 

These were comparable at high altitudes (250-400 km) but the wind affects 
dominated the 'temperature induced variations lower down* Since the two 
mechanisms are basically in opposition* it was expected that O would be 
depleted at low altitudes (<200 ka) due to the vertical wind and slightly 
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enhanced at high altitudes (>400 kra) because of thermal expansion* A 
quite different situation was predicted for Nj* Being the major species 
and diffusing negligibly through 0, the gas was found to be controlled 
mainly by thermal sources throughout the upper atmosphere although the 
circulation system was responsible for enhancing the concentration at high 
latitudes while depressing It near the equator* 

An Increase of the mass density in the thermosphere was linked to 
enhanced auroral electrojet activity by Forbes and Marcos (1973) who 
examined accelerometer data taken by the Air Force 0VL-15 satellite and 
orbital decay information on the 1968-64A Satellite* The density 
increases found in the days id o aurora agreed well with Cole's ( 1 962* 1966, 
1971) Joule heating theory* A wavelike structure originating in the 
auroral zone was also observed during periods of electrojet activity* 

Several equatorial variations of the composition were discussed by 
Reber et al* (1973) using Ogo 6 mass spectrometer data* Temperature 
enhancements of 50**130 K at low latitudes and up to 1000 K at midlatitudes 
were Inferred from M 2 density measurements taken during periods of high 
magnetic activity* Substantial decreases of n(0) and n(He) and increases 
of the N 2 concentration coincided with these magnetically disturbed 
periods* Variations of with latitude near the equator appeared to 
correspond with F-r.sgion electron density changes which had bean 
associated with the geomagnetic anomaly* The diurnal maximum of He at 450 
km was observed near 1000 LT, while N 2 and O peaked nearer to 1500 IT. 
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Haya at al. (1973) Investigated the effects of auroral heating jby 
particle precipitation and by electric fields on the thermospheric 
composition* In particular, changes of composition resulting from 
vertical winds which were excited by the heat generated by magnetic storms 
were estimated by solving the vertical continuity and momentum equations 
for the period of the sudden heating* Horizontal variations of 
composition were assumed to be negligible during the onset of a magnetic 
storm* The mean molecular mass was shown to increase at 500 km regardless 
3 of the type > s) f energy input and the largest variations of the composition 
associated with the particle fluxes resulted from 'soft' (< 0.6 keV) 
precipitation* This was explained by the fact that winds cause the 
greatest changes where the vertical composition gradients are largest* 

The region over tdiich this occurs corresponds to an area of large electric 
fields and is also the zone of 'soft' particle precipitation (near 130 
km)* Estimates of the high altitude variation of composition for similar 
heating profiles showed that this variation was nearly independent of the 

it-, 

did vary significantly with the total energy 
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injected into the system. 

The response of the thermosphere to geomagnetic heating was examined 

in detail by Mayr and Volland (1973b). The theoretical model of Mayr and 

\\ 

Volland (1972b) with a modified height distribution of the heat input was 

l 

used to estimate the variations in thermospheric temperature and density 
due to Joule heating* The N, density increased at all altitudes for every 
harmonic as would be expected* However, the net increase due to the third 
harmonic was much smaller than that caused by the first, indicating 


rate of heating although it 
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substantial damping by advectlon* The 0 density decreased at low 
altitudes (<2QQ km) as a result of upward diffusion through Nj# The 

o •( 

diffusion was caused by the higher harmonics of the heat Input (P^ through 
I» 8 ) and persisted to high altitudes, becoming offset by thermal expansion 
only above roughly 400 km. This explained the relatively small variations 
seen In n(O) at high altitudes (400 km)* The situation was completely 
different for helium* In this case, the diffusive effects controlled the 
density distribution throughout the thermosphere because of the light mass 

end low concentretlon of He* Thermal effects only slightly offset the 

6 

diffusive Influ e nces on He above 300 km* The neutral temperature changes 
were also quite marked. The first harmonic of the heating source caused 

V 

sn overall increase in the temperature but with little latitudinal 
dependence* The higher harmonica, through vertical diffusion and 

o 

horizontal transport, showed large departures from the * diffusive, 


equilibrium' temperature determined by the first harmonic. The so 
departures were strongly dspandsnt on lscitude showing the importance of 
including diffusion in calculations dealing with the global temperature 
distribution* 
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Trinks etal. (1975) Inferred from Esro 4 gas analyzer and 
ground-baaed iohoaOnde measurements west of 180° E longitude that 
midlatitude 0 to Nj ratio between 250 and 400 km substantially decreased 
during several magnetic disturbances in February, 1973* This ratio was 
fOund to decrease even more at high latitudes In a separate apd apparently 
unrelated region* The raid-latitude behavior was attributed to energy 
deposited below 150 km and transported from the high latitude auroral 
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oval*! possibly by wind* or gravity vavaa* 
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Empirical avid one a of compositions! changes during a period of 
increasing magnetic activity was obtained by Prolsa and Frick* (1976) from 
Earo 4 tnaaa spactromatar samples* The disturbed to quiet time density 
ratio R(n) served as an indicator of compositional changes throughout the 
period studied* R(n) increased above 1 for N, and for At at high 
latitudes and became especially large near the peak of the magnetic 
disturbance} R(n) decreased below 1 for He In the same region becoming as 
small as 0* 1* The behavior of 0 was more complex and both increases and 
decreases of R(n) were observed* Generally the decreases appeared at 
lower altitudes* An upper limit of one orbital period (ninety minutes) 
was placed on the high latitude response time by the authors based on the 
results of this study. 

Mayr and Volland (1976) again uaod Joule heating as an auroral energy 
input to determine what role density waves excited by the aurora played in 
modifying the composition in the thermosphere* Diffusion was found to be 
responalbl* Cor tive # large phase difference of about 220°batwaan lie and 
waves and of their relatively equal amplitude, Hie phase difference was in 
agreement with Reber at a!» (1975) although there remained some question 
as to the actual amplitudes* Th« wavelengths considered were of the order 
of 1000 km and other modes, if present, were said to be able to 
significantly altar tbaae results* 
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Neutral N 2 * 0, Ar, and He data taken down to 160 km by Atmosphere 
Explorer C during a magnetic storm in February 1974 were analyzed by Hedln 
at alt (1977c) in order to determine the response of the composition at 
low altitudes to rasgnetlc disturbances* The N 2 and Ar concentrations were 
found to increase and the 0 and He densities to decrease above 45°N 
geographic latitude* while below this latitude all densities increased 
during the storm* A very slow return to normal conditions after the 
disturbance! indicated that the atmosphere has a tlma integration effect 
with respect to this type of energy input* There was no apparent 
correlation of the neutrals to the 3-hour ap index* but a correlation with 
magnetic latitude was aeon* The authors noted that this correlation was 
asymmetrical suggesting a possible magnetic local time* universal time* of 
longitude influence* The 0 to N 2 ratio was also found to be Closely 
correlated with in-si tu electron densities and F 2 peak densities as 
measured by lonosondes* 

1. 2*4*3 Neutral Wind Observations 

Theoretical models of the global circulation system do not yet 
accurately represent the small scale features* This Inaccuracy stems from 
the global empirical pressure and temperature models used to determine the 
circulation patterns* The global model of Jacchia (1971)* which is based 
on satellite drag data* exhibits a ^lurnal pressure maximum near 1400 LT 
while incoherent scatter results (Nisbet, 1967) show the temperature 
maximum to be near 1600 LT* Studies of the Ogo 6 alrglow data (Blaraont et 
al*, 1974; Thuillieret al** 1977; Nisbet et al«* 1977) have revealed that 


the maximum temperatures occur near the local summer pole at solstice and 
not near the equator as was previously assumed. Further, Mayr and Volland 
(1972b) have determined that the temperature and mass density do not 
maximize at the same local time because of transport effects; the 
temperature appears to peak near 1600 LT while the diurnal maximum of the 
mass density is nearer to U00 LT. 


While more recent models of the thermosphere (Hedin et al., 1977a, 

b; von Zahn et al>, 1977) have improved the pressure and temperature 

0 

estimates, the wind field is still too variable to be accurately 
^represented on the basis on global models. For example, meridional winds 
of substantial magnitude have been 3hown to flow from the p^lar cap. 

Nisbet and Glendr (1977) have found that these winds are nofe uniform but 
depend strongly on magnetic local time. It is precisely this type of 
dependence that the models have yet to incorporate. For this reason, only 
empirical results will be considered here. 


% 

Kent (1970) has discussed two techniques for determining neutral 
winds. One is based on ion drift calculations, which are used to infer 
neutral velocities and the other is direct measurement of the neutral 
motion. Ion velocities can be found using Incoherent scatter facilities 
while rocket vapor trail measurements are often employed as a direct 
approach. Other direct methods have been developed and will be discussed 
later* 

0 ' 
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Reddy (1974) observed a meridional circulation pattern using 
Incoherent scatter data from St. Santln (44.7°tl) during several days in 
February* 1969, when a magnetic storm was in progress. Equatorward winds 
were seen, above roughly 120 km, while a poleward flow appeared nearer to 
100 km. Estimates of the momentum balance showed that the 100 km 'return 
flow' could support the higher altitude equatorward flux. An Interesting 
observation of the wind structure made was that -the difference in velocity 
between quiet and distrubed days peaked at 165 km and decreased above this 
altitude. The implications of this observation towards understanding the 
nature of the wind source were not known. 

Support to the circulation theory of magnetic storms was lent by 
Brekke et al. (1974) who Inferred neutral wind speeds from Incoherent 
scatter data taken at Chatanika, Alaska (ILAT-65 0 ) during the 
magnetospherlc storm of 3 August to 9 August, 1972. The magnitude of the 
daytime antisolar flow reported by Brekke et al. (1973) decreased during 
the storm and, at one point, even reversed direction. This implied that 
strong local forces at high latitudes were effectively competing with the 

r o 

predominant day-to-night flow caused by solar heating. Since ion drag was 
not expected to cause such a change, a heat source north of Chatanika 
which produced a local pressure gradient was postulated. The nighttime 
equatorward flow resulting primarily from the solar tide was reduced, 
presumably due to an increase in ion drag, while the afternoon and early 
morning winds experienced a shift in direction towards the equator during 
the storm. These effects were thought to influence the motion far from 
the auroral- zones, and it was suggested that they be included in future 
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Investigations of the global wind system* 

Eoble et al* (1977) used incoherent scatter measurements of the 
electron and ion temperatures* electron density and vertical plasma drift 
to determine neutral winds over Millstone Hill (42*6°M, 7l*5°W) for three 
winter and three summer georaagnetically quiet days in 1969 and 1970* The 
calculations were made following Roble et al* (1974)* These authors had 
deduced neutral winds over Millstone Hill for two quiet days near equinox 
and found diurnal averages of the zonal wind to be westward at about 30 m 
s“ l . In the more recent work* the nighttime equatorward winds found were 
larger for summer that for winter, but the daytime poleward winds were 
smaller during summer than in winter* The diurnal average summer winds 
were 15 m s"*^ westward for the zonal and 45 m s~* equatorward for the 
meridional component* The winter averages were 15 m s** 1 eastward for the 
zonal but a highly irregular meridional component was found. All 
velocities were calculated in a form consistent with the Ogo 6 model 
density values. These values were used in the ion heating and cooling 
rate estimates which yielded neutral temperatures based on a locally 
balanced ion energy equation* The contributions of electric fields to the 
wind system were not Included since fields on the order of 1-2 mV m * at 
this location (Kirchoff and Carpenter* 1975) were felt to cause only small 
Velocity changes. 

Midlatitude winds do not appear to depend critically upon electric 
fields, but high latitude motions are much more tightly coupled to the a 
E-fields, particularly in the auroral zones* through ion drag and Joule 
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heating (Fedder and Banks , 1972)* Therefor*, high latitude analyses oust 
take the Ion drag force Into account* 

Simultaneous measurements of Ion drifts and neutral flows at high 
magnetic latitudes (.>65°) and high altitudes (>200 km) were carried out by 
Meriwether et al* (1973) using barium clouds from rocket releases to 
monitor the motions* Optical triangulation allowed a determination of 
both drift and wind velocities In the auroral zones, polar cap, and the 
transition region* Ion drag forces were sufficient to explain most of the 
wind measurements In the evening and midnight regions while other forces, 
Including a large Inertial component, were needed to explain the morning 
flows* Spatial rotation of the wind direction from the auroral zones to 
the polar cap In the evening and midnight regions was quite noticeable* 

The westward to eastward Ion convection reversal near the Harang 
discontinuity before midnight was particularly outstanding*. 

Hays and Roble (1971) obtained meridional wind speeds representative 
of the 250 km region at Michigan (42. 3°H, 83.7°W) during two geomagnetic 

o 

storms by measuring the Doppler shift of the 630 nm line of atomic oxygen. 
Southward velocities of 250 to 400 ras”* were observed throughout the night 
during both storms. A combination of diurnal and enhanced auroral zone 

sources were suspected to be responsible for these Incense flows. 

0 .... 

Nagy et al. (1974) also measured the Doppler shift of the 630 nm line 
of atomic oxygen from Ester Done, Alaska (IlAT-65 0 ) In order to Infer 
neutral wind velocities while simultaneously measuring Ion drifts by 
Incoherent scatter radar at Chatanika, Alaska. The Ion drift reversal 


near magnetic midnight seen by Meriwether et al* (1973) was quite 
prominent and the westward Ion drift In the premidnight region appeared to 
control the neutral zonal flow* although no such control lias apparent In 
the neutral meridional winds since the north-south Ion flux Is much 
smaller* The postraid night magnetic local time region displayed virtually 
no neutral zonal motion and there appeared to be no simple explanation for 
this* The authors noted In conclusion cipb a thorough Investigation of 
the wind system would require estimates of several driving 
forces: electric fields of magnetospherlc origin* diurnal pressure 
gradients from solar heating* Joule and particle heating in the auroral 
zones* and Coriolis terms* 

Monthly variations In nighttime neutral winds and temperature near 
250 km during solar minimum were catalogued by Hernandez and Roble (1977) 
using airglow data taken at Fritz Peak Observatory (39.9°n, 105. 5°W) from 
April* 1975 through March, 1976. Measurements during quiet geomagnetic 
conditions were compared to values calculated from Ogo 6 model pressure 
gradients (lledin et al., 1974). The zonal flow was predominantly eastward 
during winter nights while shifting westward during the mornings. This 
transition became progressively earlier towards summer and approached 
midnight at the solstice. The diurnally averaged flow was eastward (35-50 
n s“ l ) during winter and westward (25-30 m s“*) in summer. The equatorial 
meridional winds were largest near midnight in summer and smallest during 
winter. The daytime poleward flows maximized in winter and minimized in 
summer. Diurnal averages of these were 40 m s" 1 equatorward and 25-35 m 
s" 1 poleward during summer o and winter* respectively. The nighttime 


38 


temperatures were lowest In January and highaat during May, July, and 

August $ Agreement between measured and calculated wind speeds was 

\) 

generally good; In fact, the diurnal averages of the calculated wind 
speeds were taken aa representative of the measured diurnal averages* 
Nighttime temperatures did not agree so well* The measured values were 

^ . If 

100-130 K higher than model estimates for summer and equinox and 50-75 K 

higher for winter* Daytime differences were unknown* The authors 

" ’ //' 

suggested that di/.ferences between measured and calculated motions might 
be reduced by Including an auroral heating term In the wind system model, 
the magnitude of this term to be determined as that which would bring 
calculated and measured values Into agreement. 

An Interesting new way of measuring neutral wind velocities based on 
mass spectrometer sampling was developed by Knutson et al* (1977). 
Vertical velocity estimates were obtained by using the open source 
spectrometer on the AE-C satellite* A determination of the periods 
between successive maxima of the spectrometer signal, which were quite 
sharp, as the satellite spun allowed the wind component perpendicular to 
the spin plane to be approximated* Horizontal wind contributions were 
eliminated by using only low altitude measurements since In this case, 
only vertical motions were perpendicular to the spin plane* The largest 
vertical wind variations were found in the post magnetic midnight auroral 
region* Velocity changes as large as 65 m s*~* were observed although no 
simple relationships with altitude or magnetic activity were obtained due 
to insufficient data* These values were not always consistent with 
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tracking measurements, but the tracking data were restricted to dawn and 
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dusk* The authors no tad that their method did not suffer these 
* & 
restrictions and that horizontal winds could be found by placing the 

satellite spin axis in the vertical direction* 
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CHAPTER II 

DATABASE OP THE MODEL AND COORDINATE SYSTEM 


2* 1 Description of the Satellite and OperatlnR Conditions 

The Atmosphere Explorer-C satellite was launched into a 68.1° 

' if 

Inclination orbit on December 16, 1973 with an Initial apogee of 4300 km 
and an initial perigee of approximately 157 km* The initial perigee 
latitude changed by .5° per day from 10°tl and the local time of perigee 
decreased by about 8 minutes per day from 1430 hours local time due to 
orbital precession. The satellite was maintained in this highly 
elliptical orbit until November, 1974, when the orbit was made circular at 
altitudes between 220 km and 320 km* The orbital period was then about 90 
minutes so that successive passes were separated in longitude by about 22 
degrees* Until this maneuver, numerous low altitude excursions were made 
in order that data could be taken to altitudes below 140 km* 


The AE-C neutral density measurements used in this study were 
obtained from January, 1974 to June»-1975, a period which corresponded to 
the low activity portion of the ll-year solar cycle. The daily burst 
adjusted 10*7 cm solar radio flux at 1 A.U* was between 66*5 and 144*7 x 
10" 22 W a" 2 Hz”* throughout this period* 
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2. 2 Instrumentation - The Ooan Source Mass Spectrometer 

* ' ;* 

Several neutral mass spectrometers were included In the payload of 
AC-C* Only the Open Source Spectrometer (OSS) vac used In the praaent 
atudy in order to maintain consistency throughout the database* This 
instrument has been described in detail by Nler at al. (1973)* Briefly, 
the spatial resolution of the OSS vna 3.2 km along the orbital path for a 
sampling rate of 16 samples per second, the normal operating rate* The 
accuracy of the number density measurements was estimated to be within 
+20% and much better relative values of the concentration were presumed* 
Gases naaaured by the Instrument included N,» 0, 0^, He, Ar and several 
minor species* 

2. 3 Data Selection and Corrections 

The OSS instrument generated several thousand spectra for each gas 
aampled during one 30-rainute perigee pass of the satellite in the 

elliptical phase* This?/ 1 amount of high resolution data was considerably 

/ 

larger than desired for the present purposes; therefore, a way of unifying 
the measurements was needed* The original database of the AE satellite 
series included Telemetry (TM) files, Geophysical Unit (GU) files, and 
Unified Abstract (UA) files. The TM files consisted of satellite orbital 
information and raw Instrument voltages. The GU files contained high 
resolution data, which consisted of spectra in the case of the OSS. The 
UA files were made up of 15-sccond averages of the high resolution spectra 
'x^d as such were well suited for the present study. By definition, these 
averages were centered on the universal time specified for the measurement 
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so that th« sotsllits Was vrithin about; 65 ktn (along the orbital path) of 
the position corresponding to this universal time when the samples 
comprising the averages were taken . 

Measurements of tig and He by the OSS were accomplished in a fairly 
straightforward manner since both of these gases ere reasonably inert* 
Estimation of the 0 density was much more difficult due to the 
interactions that this highly reactive gas underwent on the surfaces of 
the instrument* In fact, in the normal mode of operation, 0 was actually 
measured as it would have been In a closed source instrument* That is, as 
0,,, after recombination on the walls of the OSS ionizing region* 
Unfortunately, the recombination of 0 into Og was incomplete and the 
measurements had to be scaled to reflect this* A correction for ambient 
Og In the atmosphere was also made* Both corrections were applied as 
described by Kayser (1977). More information concerning the analysis and 
calibration of the OSS as well as techniques used to measure 0 and Og via 
mass spectrometers have been supplied by Nier et al. (1974, 1976), French 
et al* (1975), and Kayser and Potter (1976). 

Lower limits were set on the concentrations as suggested by Kayser 

(1977) to reduce the problems of outgassing and low multiplier counting 

12 -3 11 

rates at high altitudes. The minimum values used were 5 x 10 m >10 
—3 10-3 

m , and 5 x 10 m for 0, He, and Ng, respectively. The altitudes at 
which these concentrations were reached depended strongly upon the time of 
sampling, the geographical position of the satellite, the magnetic and 
solar activities at the time the measurements were taken, and the neutral 


gas sampled* 

Nb o 
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2.4' The Marcnetlc Coordinate System 

; .. 7 

Many aeronomic studies have organized data in coordinate systems 

fj 

v»hich are based on she magnetic field of the earth* The 'geomagnetic* 
coordinate system is derived from a centered dipole approximation to the 
magnetic field./, This representation is equivalent to the superposition of 
two dipoles located at the center of the earth, one aligned along the 
geographic axis and the other lying in the equatorial plane. Geomagnetic 
latitude is defined analogously to geographic latitude* The geomagnetic 
equator is defined as the great circle around a spherical earth which is 
normal to the centered dipole. Geomagnetic local time is defined in terms 
of the angle formed between the plane containing the earth-sun line and a 
vector parallel to the dipole, and the plane containing a vector from the 
center of the dipole to an observer and a vector parallel cro the dipole. 

'Magnetic', 'eccentric dipole', or 'offset dipole' coordinates are 
derived by displacing the centered dipole some ^distance from the center of 
the earth in order to minimize the second term in the harmonic expansion 
of the magnetic field without affecting the first term. Magnetic latitude 

, U • 

and magnetic local time are then def ined^ analogously to their geomagnetic' 3 
counterparts. 

The transformation from geographic to geomagnetic coordinates is 
simpler than that from; geographic to magnetic or eccentric dipole 
coordinates. For this reason, geomagnetic coordinates were preferable to 
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magnetic coordinates before computers were widely used* Eccentric dipole 

f \ 

coordinates^ however, are much more closely related to the magnetic field 

W " o 

configuration at high latitudes, in the auroral zones, and in the polar 
cap* Because they yield a better representation, eccentric dipole 
latitude and magnetic local time will be used exclusively in the present 
analysis • 


The equations used to transfer from geographic to magnetic 
coordinates were given by Agy (1965) and appear in Appendix 1 of 
Matsushita and Campbell (1967). These equations used the direction of the 
centered dipole derived by Finch and Lenton (1957) and the displacement of 
the eccentric dipole calculated by Parkinson and Cleary (1958) • Secular 
variations in the magnetic field cause quantities to change with time and 

It was found desirable to update them to more recent values* The I.G.R.F. 

0 

1975 magnetic field model in conjunction with the method described by 
Bartels (1936) was used to do this. The updated equations for magnetic 


latitude and magnetic local time nrf 

Magnetic latitude ■ sin (§ • ft) 

Magnetic local time ■ (tan + 180°) 


( 2 * 1 ) 

( 2 * 2 ) 


wheje 

S’ ■ S x ft • vector in magnetic equatorial plane perpendicular 
to earth-sun line * ii 

2T - | x ft * vector in magnetic equatorial plane perpendicular 
to magnetic position vector 
ft - unit vector in direction of (negative) dipole 
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- i cos 289.5° cos 78.6° + 5 sin 289.5° cos 78.6° + 

*» 

ft sin 78.6° [based on I.G.R.F 1975 J 
f d cos 146.9° cos 19.5° + J d sin 146.9° cos 19.5° + 

£ d sin 19.5° [baaed on I+G.R.F. 1975] 
d ■ magnitude of displacement: of eccentric dipole 
» 474 km [baaed on I.G.R.F. 1975] 

§*■ position vector relative to dipole (magnetic position 
vector) 

■ |T - <T 

pr - position radius vector 

• it cos X cos <|> Hh 5 r sin X cos $ + & r sin <|> 
s - unit vector in direction of (infinitely distant) sun 

- i cos (180° - uT) cos 6 + J ain (180° - wT) cos 6 + 
k ain 6 

£,$,£■ unit vectors from center of spherical earth along 
direction indicated! i to intersection of equator 
and Greenwich meridian, j to intersection of 
equator and 90° east geographic 
meridian, k to north pole 
R ■ radius of earth - 6371.2 km 
h * height above the earth, km 
r - R + h 

<j> - geographic latitude 
X • geographic longitude (east) 

T ■ Greenwich local solar time, hours 


■- *V 

■U', 


0 


fp 


Vi'. , 

P 


***.*■■34 'M 






w • 15 d«8 h"** (180° - ujT --Greenwich hour angle of sun) 

6 - solar declination 

_ In these equations, the position used for the geomagnetic north pole 
was 78«6°N, 289* 5°E and the displacement of the dipole occentrlc from the 
center of the eerth was calculated to be A 74 km towards 19«5°H, 146»9°E. 
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CHAPTER III 


EXAMINATION OF THE LONG-TERM TRENDS OF THE COMPOSITION IN THE MAGNETIC 
COORDINATE SYSTEM 

I 

3* 1 Determination of Average Density Structure 

The behavior of the neutral thermospheric composition has sometimes 

o 

been described on the basis of data taken from individual satellite orbits 
or groups of orbits* Such limited Information must be interpreted very 
carefully, however, since temporal variations due to magnetic storms and 
large gravity waves can cause large orbit to orbit variations which mask 
the global and seasonal variations* 

One way to unmask the basic structure is to reduce these local 
effects by analyzing data taken over long periods of time which have been 
corrected for long-term solar and geophysical variations* These data 
should then be more representative of average conditions making it much 
easier to separate persistent from transient distributions of the neutral 
gas more clearly. In order to study the density structure in the magnetic 
coordinate system and to find persistent effects this approach will be 
followed here* 

Several solar and geophysical influences on the composition need to 

be considered before averages are taken* The contribution of solar 

activity variations towards neutral density changes during 1974 and the 

first 6 months of 1975 did not seem to be of great importance. The 

11-year solar cycle was near minimum and the F 1rt , 

lu* 7 


index generally changed 
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by lean than 40 x 10 W ra Hjs • Baaed on the MSIS model (Hedin et 
al», 1977a, b), It can be estimated that auch changes in solar flux 
translate into less than 10% variations in neutral temperature and 

o 

density* These fluctuations are much smaller than those arising from 
other major sources and were,°therefore, neglected* 

The influence of the level of magnetic activity on the neutral 

atmosphere, on the other hand, was found to be extremely important and was 

included in this study* The 3-hour planetary ap index was used as an 

indicator of magnetic activity because it applies globally and is readily 

available, although it is by no means the best indicator of auroral 

disturbances as will be shown later* Three levels of magnetic activity 

were considered and delineated by the following ranges of ap values: 0 to 

4, 5 tta 15, and 18 to 67* These ranges were chosen to be representative 

of low, medium and high magnetic activity conditions, respectively. It 

should be noted that there were fewer data and more variable data taken 

over the high range of ap but much meaningful information is still 

obtained for these conditions* The dependencies of the atomic oxygen and 

helium data on altitude were treated individually and Che reasons for this 

are described below* In each case, all data were reduced to a common 

. .0 •. 
altitude (Sections 3.2, 3*3) to aid in representation and interpretation. 

Later, these fixed altitude results are generalized* 

The period over which all measurements were treated as a single 
(data) set was quite long.” All neutral density values taken within 45 
days of a given equinox or solstice were chosen to comprise a 'season'* 
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In order to determine the appropriate corrections for annual effects in 
the densities taken at the extremes of these periods, scatter plots were 
made of the data at fixed altitudes within a small range of magnetic ,, 
latitudes (+5°) and magnetic local times (+3 hours) as a function of day 
number for the entire database* 

» .. 

Tire helium data showed a winter pole to summer pole ratio of about AO 
to 1 at 300 km (Figure 1). This large variation was also quite regular 
permitting a logarithmic correction to the equinox data to be made based 
on the day number of the measurements and the magnetic latitude band into 
which the data fell. The solstice helium data were not corrected in this 
way since a deviation from uniform coverage for these seasons would have a 
much smaller effect than at equinox* The scatter in the atomic oxygen 

measurements was greater than any annual variations which might be present 

0 

(Figure 2); therefore, no corrections were made to these data. 

Finally, the magnetic coordinates at which the measurements were 
taken were considered by forming 5° 'bins' of magnetic latitude and one 
hour bins of magnetic local time centered at 5° intervals from the equator 
and on the hour, respectively. All data taken within a particular range 
of ap and corresponding to a given season falling into the same bln were 
averaged. The number of points in each bln ranged from 0 to over 100. 

Bins which had no data, including those poleward of about 78° MLAT in the 
northern hemisphere and 82° MLAT in the southern hemisphere, where the 
satellite did not reach, were filled by logarithnic interpolation. Once 
this global 'grid' was established, contour plots of the data in each 


Summer Fall Winter Sprin 



Figure 1 Scatter plot of logjQ n(Ue) at 300 km as a function of day 

number. All tfata^were taken from 50 - 60° magnetic latitude 
and between 1200 and 1800 hours magnetic local time. Similar 
plots for other latitude ranges were lused to derive a correc- 
tion factor for helium data depending on number of days from 
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hemisphere ware made. Thesa plots ata displayed in polar Corn as 
functions o£ magnetic latitude and magnetic local tima for various seasons 
and ap rangas and will ba described in detail in Section 3* 4* 

3* 2 Reduction of Atonic Oxygen to 120 km 

Atonic oxygen is the main constituent in the middle thermosphere and 
at altitudes of the order of 300 km winds cause the horizontal density 
gradients to ba small* This results in vertical transport in regions of 
energy input where the exospheric temperatures are high* In order to 
separate the effects of scale height variations due to temperature 
variations from those due to vertical diffusion at low altitudes it is 
convenient to use an equivalent atomic oxygen density at 120 km that would 
reproduce the observed variations at satellite altitudes assuming a 
diffusive equilibrium profile* these densities are a sensitive indicator 
of Vertical transport and heating in the thermosphere and provide simpler 
boundary conditions for a model* Care must, of course, be taken not to 

S .5 

interpret these literally as measurements of the actual atomic oxygen 
densities at 120 km which nay be different due to the effects of relative 

, J / 

diffusion above 120 km as discussed by Nisbet and Glenar (1977)* 

These authors have shown that measurements of the atonic oxygen and 
molecular nitrogen density at any altitude may be combined with a neutral 
temperature in a straightforward manner so as to give estimates of the 0 
density at 120 km if an N, density and a neutral temperature are assumed 
at this lower boundary* They have also shown that the resulting 0 values 
at 120 km are relatively insensitive to the value used for the neutral 


temperature at the height at which the measuremente were taken (leas than 
a AS change in density Eor a 100 K temperature change)* This result 
allows model temperatures to be used in place of actual measurements with 
very little change *■’ accuracy* Their method will be followed heret 

\f' 

The density of a neutral gas in diffusive equilibrium at any height * 
above a certain base level z may be found ignoring thermal diffusion 

Q ° 

effects, which are negligible for 0, from the equation, 

T(*J , , «t,\ 

"i<*> ■ "!<*„) tfnr 1 * x p< '-r> / Z(J d,> (3,1) 


where 

z level where density is desired, km 
llj molecular weight, g mole”* 

R gas constant, J mole” 1 !?” 1. 


Using in c hi» equation, the Integral can be written as, 



<2> 


Substituting this into the first equation with 0 instead of as the 

gas. 











1 
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ii [L, V*/ * 

In n 0 (z o ) - In n 0 <z) - (-.S-j i„ [-~^-yl + (g -2. .[] ln t“^|j3 (3.3) 


Now, for a base level z o - 120 km and N Q - 16, ^ - 28, this 


equation becomes,. 


1<>8 i0 V 120 > " lo «io n o (l) + 7 Ilog lO n l|, < l20 > 

- lo S l0 n, (*)] + | lo Bl0 


(3.4) 


' ^ ° ' %l * 0 ' 

where the cqbbqii. logarithm of the density has bc«f» mm . 

,, . • & 

S'- ( 

If Che Mj density and neutral temperature are taken from die CIRA 
(1972) 1000 K model <n„ -3,79 a 10 17 m’ 3 . T -334.5 K). . * 

los 10 n O (l20 > " lo 2 io n 0 (z> + 7(17.579 - log^^ (*jj 
+ | log.^ rt.i5L.i <, 

7 iUS 10 l 334.5 J (3.5) 

0 

0 , v A . 

„ .. s' 7, 

/( a u ci ' ■' ,v •• . . .9 , 

Thia equation was used, to reduce the 0 data to 120 kra. 

3»3 Reduction of Helium to 300 km <. 

Helium Is one of the lightest constituents In the thermosphere and as 
such has a very large scale height. Under the assumption of hydrostatic 
equilibrium, the global helium distribution will be fairly similar at most 
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I,:' 

| 

• 

k 

I 



thermospheric heights* Therefore, an intermediate altitude such aa 300 km 
serve a as a convenient reference height to study the variations in neutral 
helium* Further* aince satellite measurements need only be extrapolated 
over leas than about one scale height to reach this altitude* it may be 
expected that the inferred distribution at 300 km will reflect the 
distribution over the height range of the satellite measurements fairly 
veil* At this altitude* the thermospheric temperatures have almost 
reached exospheric values* The assumption of an isothermal atmosphere is* 
therefore* reasonable and will be made* 

The neutral density of He in diffusive equilibrium Is given by 


, . / \ r T{z) , — t t “**He f s s(z)dz 

n l!e " "lle^o 5 { ~t7z“) ®* Pt J z q "f?z7" 


f Z dT y 

" “/ Z o TU ) 1 

where a is the thermal diffusion factor, 
-•38 (Banks and Kockarts* 1973) 


(3*6) 


In an isothermal atmosphere* this reduces to 

”^He z 

V (8) " n ne (t o > «P< Hr~ ] f z S(*)da> 0-7) 

w 0 


where T w is the exospheric temperature. 
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After integration and substitution of z 0 * 300 km into the above 
expression, the helium density at 300 km nay be calculated as, 

H Hb SC°> 

n |Je (300) - n„ e (z) exp^^— (b(a) - h(300)J> (3*8) 

where 

h(z) la 3 eopotentinl height of z, kra 

g(0) is acceleration of gravity at earth's surface, 

9.81 n sec - * 

This expression and the JISIS model neutral temperature (Kedin at al. , 
1977a) corresponding to the position of the satellite and the solar and 
geophysical conditions present for a given measurement vers used to 
determine the 300 km helium densities. 

3.4 Comparison of the Density Averages Under High and Low Solar 

Activity Conditions 

3.4*1 Atonic Oxysen Densities at 120 km 

Long-term averages of the atonic oxygen densities for spring equinox 
conditions are shown in Figures 3, A, and 5. Data from the AE-C satellite 
were used to moke the low solar activity contours (top) and Ogo 6 
satellite measurements (Stehle et al*, 1979), which represented high solar 
activities (Fj^ y ■ 101 to 239), were used to construct the contours for 
the solar maximum (bottom) . The density variations for low magnetic 
activity levels, ap - 0 to A, are illustrated in Figure 3. Atonic oxygen 
shows similar high latitude structure for both high and low solar 
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activities • The minimum of n(0) at 120 km extends from the magnetic po^es 
through the postmidnight magnetic local time (MLT) sectors in both cases, 
although this minimum is less restricted for Ogo 6. This is probably due 

O 

to fewer measurements in this region by the Ogo 6 satellite. The 
densities in these minima axe r abbut^50% of the respective hemispherical 
averages, but the absolute values of tl^e AE-C data are nearly 40% lower 
than the Ogo 6 measurements. This is in\agreement with Hedln et 
al. (1977b), who predict an increase in n(0) with increasing solar 


activity. The high latitude structure becomes sharper for medium magnetic 

activity levels, ap ■ 5 to 15, and the density minima persist in the 

<2 

postmidnight MLT region under both high and low solar activities (figure 

o : ■ # 

4) s> The densities in these minima, become smaller than for low magnetic w 

activities, being only about 40% of the hemispherical averages. The 

atomic oxygen densities show even more structure for high magnetic 

activity levels, ap * 18 to 67 (Figure 5). For both solar activities, the 

densities are smallest near 0700 MLT in the polar regions. The values in 

these minima are less than 30% of the hemispherical' averages and only 

j . ,, ; v 

about 50% of the corresponding low magnetic activity values. The 
postmidnight MLT density minimum persists for the summer solstice atomic 
oxygen dath (Figure 6). For this season, the AE-C minimum densities are 

about 30% of the northern hemispherical average but the Ogo 6 values are 

© J 

nearer to 40% of this average, suggesting a stronger annual effect for low 

r j 

solar activity condition. The Winter n(Q) values (Figure 7), while 

showing somewhat different MLT structure near the magnetic poles, support 

./■ ' /(' ....... ...... 

this a3 the winter to summer pole-to*?pole ratio for AE-C is about 2 to 1 
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Acomic oxygen density at 120 km for spring equinox 
and ap values between 0 and 4 as a function of mag- 
netic latitude and magnetic local time. Ton repre- 
sents low solar activity conditions (1974) and AE-C 
data and bottom shows high solar activity conditions 
(1970) and Ogo 6 data (Stehle et al., 1979). 
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Figure 4 Atomic oxygen density at 120 km for spring equinox 
and ap values between 5 and 15 as a function of 
magnetic latitude and magnetic local time. Top 
represents low solar activity conditions (1974) 
and AE-C data and bottom shows high solar activity 
conditions (1970) ar.d Ogo 6 data (Stehle et al., 
1979). 
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Figure 
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Atomic oxygen density at 120 km for spring equinox 
and ap values oetween 18 and 67 as a function of 
magnetic latitude and magnetic local time. Top 
represents low solar activity conditions (1974) 
and AE-C data and bottom shows high solar activity 
conditions (1970) and Ogo 6 data (Stehle et al., 
1979). 
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Figure 6 Atomic oxygen density at 120 km for summer solstice 
and ap values between 5 and 15 as a function of 
magnetic latitude and magnetic local time. Top 
represents low solar activity conditions (1974) 
and AE-C data and bottom shows high solar activity 
conditions (1970) and Ogo 6 data (Stehle et al., 
1979). 
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Figure 7 
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Acomlc oxygen denslcy at 120 km for winter solstice 
and ap values between S and 15 as a function of 
magnetic latitude and magnetic local time. Top 
represents low solar activity conditions (1974) 
and AE-C data and bottom shows high solar activity 
conditions (1969) and Ogo 6 data (Stehle et al.. 



but nearer to 1*5 to 1 for Ogo 6 


The atonic oxygen densities decrease at high latitudes for increasing 
magnetic activities under both high and low solar activity conditions* 

They are also lower in summer than in winter for both cases and show 
similar and persistent structure in the region of the westward electrojet 
where large density variations are expected due to Joule heating from 
ionospheric currents. That this structure is maintained over one-half of 

o 

the 11-year solar cycle suggests that the magnetic coordinate system is 

O 

well suited for modeling the neutral composition* 

3*4*2 Helium Densities at 300 km 

The averages of the helium data for AE-C and Ogo 6 during spring 
equinox are represented by Figures 3, 9, and 10. The low magnetic 
activity, ap - 0 to 4, structure of helium at high latitudes is not 

o O 

apparent in these contours probably due to a paucity of data, especially 

o -- 

for Ogo 6. At medium levels of magnetic activity, ap ■ 5 to 15. the 

O 

spring densities show high latitude minima near the magnetic pole for AE-C 
and over the magnetic pole for Ogo 6. The n(He) values in these regions 
are less than 20% of their respective hemispherical averages. These 
density depressions become more pronounced at high magnetic activity 
levels, ap - 18 to 67, although they are no smaller in magnitude than the 

c> ■ . . . ' 

densities for medlumap. This 'saturation' of density variations with 
increasing ap has also been sean by lie din et al • (1977c) in AE-C data 

taken during a magnetic Strom in February, 1974. The helium 
concentrations in the postnidnight minima for high ap levels decrease to 

rt . ... 0 
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about 40% of the low magnetic activity values for both quiet and active 
solar conditions. 

Summer solstice heliun averages (Figure 11) maintain the high 
latitude equinox structure having minima in the postmidnight MLT area 
under high and low solar activities although the Ogo 6 density depression 
extends into the premidnight MLT region near the eastward electrojet. The 
winter solstice data (Figure 12), on the other hand, maximize over the 
magnetic pole and a very strong annual effect seems to be the reason for 
this. The winter to summer pole-to-pole ratio for the AE-C helium data 
is greater than 25 which is close to that reported by Mauersberger et 
el. (1976) and Straus et al. (1977), the latter from a theoretical model 
Including dynamical effects. This ratio is less than 10 for the Ogo 6 
data and an increase in exospheric transport with increasing exospheric 
temperatures and solar activity as suggested by Reber and Hays (1973) and 
Strauss et al. (1977) would explain this since smaller density gradients 
in latitude would be maintained in this case. 

The helium density averages and the atonic oxygen density averages 
show very similar behavior at high latitudes under both quiet and active 
solar conditions for nearly all levels of magnetic activity, except during 
vinter when the large annual variation of helium obscures this 
structure. These two constituents have masses less than molecular 
nitrogen and are very good Indicators of horizontal winds, 
thermospheric heating and upwelling. Since their structure is much he 
same in the magnetic coordinate system under widely varying solar and 
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Figure 8 Helium density at 300 km for spring equinox and ap 
values between 0 and 4 as a function of magnetic 
latitude and magnetic local time. Top represents 
low solar activity conditions (1974) and AE-C data 
and bottom shows high solar activity conditions 
(1970)and Ogo 6 data (Stehle et al., 1979). 
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Figure 9 


Helium density at 300 km for spring equinox and ar 
values between 5 and 15 as a function of magnetic 
latitude and magnetic local time. Top represents 
low solar activity conditions (1974) and AE-C data 

* hoWS high • olar activity conditions 
U970) and Ogo 6 data (Stehle et al.. 1979} . 










MAONCTlC local timc 


Figure 10 Helium density at 300 km for spring equinox end 
sp values between 18 and 67 as a function of mag 
netic latitude and magnetic local time. Top rep 
resents low solar activity conditions (1974) and 
AE-C data and bottom shows high solar activity 
conditions (1970) and Ogo 6 data (Stehle et al.. 
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Helium density at 300 km for sunnier solstice and 
ap values between 3 and 15 as a function of magnetic 
latitude and magnetic local time. Top represents 
low solar activity conditions (1974) and AE-C 
data and bottom shows high solar activity conditions 
(1970) and Ogo 6 data (Stable et al., 1979). 
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Helium density at 300 km ,'or winter solatlce and 
ap values between 3 and 13 as a function of mag- 
netic latitude and magnetic local time. Top 
represents low solar activity conditions (1974) 
and AE-C data and bottc' shows high solar activity 
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geophysical conditions, it would appear that a general analytic nodal of 
their distribution in this ays tarn would be of value in studies of energy 
deposition processes and dynamics of the polar thermosphere. Such A modal 
,cloped in the next section (Chapter IV), 

O 




CHAPTER IV 


METHOD OF ANALYSIS AND CONSTPvUCTION OF A MODEL 


4. 1 The Neutral Magnetic Coordinate Model 

• 

The long-term averages of the atomic oxygen and helium densities 
showed consistent structure at high latitudes for various seasons and 
magnetic activity levels. Extraction of annual variations from these 
averages is difficult due to the long time periods used and determination 
of semiannual variations, if any, is impossible. In order to include 
these effects in the representation of the atomic oxygen and helium < 

if 

distributions. Individual data points were used to construct the global 
model. Enough measurements were included (over 2500) that the large and 
highly localized temporal variations discussed earlier should not 
adversely affect th-ji results. Of course, all small scale effects will not 
be faithfully reproduced by a model of this type, but a wide range of 
solar and geophysical conditions will be covered globally. 

4.2 Temperature Distribution 

Molecular nitrogen is the major constituent in the region between 90 
km and 150 km where relative diffusion of the atmosphere constituents is 
of importance.,, The mass of N.g is very close to the mean molecular mass in 
this region which means Chat is not greatly affected by the transition 

. . i'} A 

from the realm in which the distribution is controlled by eddy diffusion 
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to the area »here free molecular dlffuaion domlnatea. The Kj altitude 
profile can, therefore, be approximated at all heights by the character- 
istics of a hydrostatic distribution even when large velocity fields 
are present. The approximation allows a neutral temperature profile to 
be inferred from the ^ measurements. In particular, each ^ density 
value was used to calculate an equivalent exospheric temperature, refer- 
red to as NT. Only data taken above 190 km were used so that a unique 
exospheric temperature could be inferred from each density value. 
Since the database of AE-C was well over 100,000 points, it was found 
Impractical to use the entire data set and only one of every ten admis- 
sible data points was actually used in calculating NT. The locations 
of these points in the magnetic coordinate system within 45 days of the 
equinoxes and solstices are shown in Figures 13 and 14. 



A Bates (1959) temperature profile was assumed for calculating NT 
and it was necessary to assume lower boundary values for the temperature 

O 

and the N, density and a shape parameter for the temperature profile in 

0 5 

order to determine NT. These values were taken from the MSIS model (Hedin 

et al . , 1977a). NT was computed by Iteration with the aid of equation 

4.3. Since the molecular nitrogen data used here were obtained at times 

K 'of low solar activity, it was difficult to estimate the effects of EUV 

variations on the temperatures. The parameters reflecting; these effects 
( ■ : f) 

were also taken from the MSIS model which va| base*/ on data taken during 

-various phases; of -'the 11-year solar activity cycle. 


o 
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Figure 13 Data used in the global model depending on magnetic latitude and magnetic 
local time. Points shown were taken within 45 days of the equinoxes and 
above 190,1cm for Ap values less than or equal to ten. Northern hemisphere 
represent^ local spring. 
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The expression used for the neutral temperature profile* taken 
largely from the MSIS model, was of the form* 

T(z,L) - T m (L) - CT„(L) - T l2Q ] exp(-s [z-12oT> 


1JL) - T G (L) 

J o 


(4.1) 


The spherical harmonic expansion, G^,(L), can be written, with brief 
descriptions of the terms involved, as 


G t m 1 
^ 0 

Time Independent Terms 


(4.2) 


+ a 10 F 10 + a 20 P 20 + a 40 P 40 
Solar Activity Terms 

+ £ 00 4F + £ 00 <4f,2+? 3o 4? 

Annual Terms ,, " . 

+ Cqq cos n(t d - t Q0 ) + (c 10 P 10 + C 30 P 30^ 

(1 + f^AF) cos n(tj - tjj) 

Semiannual terms 

+ (c 00 + C 20 P 20 ) C0s 2fJ<t d " + C 10 P 10 C0S 2n(t d *10* 

Diurnal Terms 


+ Ea ll P 1 1 + 3 3 1 P 3 1 + a 51 P 51 + ^ C 11 P 11„ + C 21 P 21^ C0S n(t d “ C 10^ 
(1 + f’jjAF) cos wt " 

(P + + b 31 ? 31 + *91*51 + (d li f ii + 4i*ll> cos a(t d - ‘to” 

(1 + f^AF) sin idt / + ]- 


o 
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Semidiurnal Terms 

4- « 22 p 22 <l + ^1 A ^ coa 2ut+ b 22 P 22 (l + ^n AP > 8in 2WT 
Magnetic Activity Terms 


c a 
00 c , 


+ + ^20 P 20 ^ A A 


where 

AA 


(Ap - 4) 


AF ■ F, a - - F.a - 
10.7 10.7 


AF 


P I0.7 * 150 


P nn associated Legendre polynomials equal to 

(1 - x 2 ) ra/2 / 2 n n! (d n + ra /dx n + m ) (x 2 - l) n , 
x ■ cosine of magnetic colatitude. 


n 


-1 


jtt - 2x7 365 day* 
id ■ 2ir/24 hour”* ■*, 
t ■ magnetic local time, hours 
L list of solar and geophysical parameters 
z altitude, km 
day of year 

daily magnetic index /( 


d 

Ap 


FiO. yl0.7 cm solar flux on previous day, 


10" 22 U m“ 2 Hz -1 


F ^average of F^ y over 3 solar 

rotations centered on t . 

e. ‘ d 

T„ average exospheric temperature, K 
T temperature at height z 
T 120* ower boundary temperature, 386 K 








78 


-2 -1 

s temperature gradient parameter, 2*548 x 10 km 


The coefficients of the spherical harmonic expansion- for T W (L) were 
determined by calculating a nonweighted least squares fit of NT to 
¥^(1,) • This fit was done in two steps. The first step included only 
data where Ap<10 and all coefficients except those due to magnetic 

Cj 

activity were determined* The second step used all data where Ap<130 (one 
third of those are above an Ap of ten) to determine the magnetic activity 

coefficients. The coefficients and their standard deviations are listed 

o 

in Table 1. " 


4* 3 Atonic Oxygen and Helium Density Distributions 


The global distribution of the neutral composition be represented 
over a range of altitudes once the temperature profile has been determined 
and lower boundary densities have been specified. The extrapolation in 

U _ 

altitude is made through the diffusive equilibrium profiles of Walker 
(1965), given by 


n(z,L) - n 120 (L) D{z,s,T 120t T (lo (L)l 

/jr X> 

O 

Q[z,s,T 120 , T^aoj - [(l-a)/(l-ae " 8? )] 1 +Y+ “ exp(-s£) (4.3) 

O 

(z - 120 ) (R e + 120) 


0 


where 
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n 120^ "^120 CC M 0-) - U 

Y - Mg 120 /sRT (L) 

8 j 20 " 9*44663 n s’" 7, 

.3 

n ambient density* m 
tfj2o overage 120 km density 

for N 2 , n u0 - 2.81 x 10 17 nf 3 and G 28 - 1 

O „i 

M mass number of constituent* 3 mole 'n 


order to determine the atomic oxygen distribution at the 120 km 
lover boundary, the n(Q) values at satellite altitudes were reduced to 120 
km using corresponding N 0 values as discussed previously (Section 3.2). 

The neutral temperatures, T * were Inferred from the N. data rather than 

Z v 

calculated from the MS IS model. The 120 km helium densities were found 
similarly and the effects of thermal diffusion for He were included in the 

reduction expression, 

r; ~ T% . 0 


In n He (120) - ln ; n He (r) + }[ln n,^(l20) - In n^(z)] 


+ tl + «J ln(a---J 

a 120 


(4.4) 


The 120 km boundary densities were used to determine the coefficients 

(Table 2) of the spherical harmonic expansion G W (L) , chosen to be 

M 0 

identical in form to that used for the inferred exospheric temperatures 


o 
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Table 1, Spherical harmonic coefficients for the inferred 
exospheric temperature 


Coefficients 

T 

« 


Coefficients 

\ . 

Average 

1 . 045E+3i3 . E+0 


4 ' 

-3.007E-2i2.E-3 

a io 

“ * 

-1.271E-2i%E-3 


c2 

c 20 

-2. 795E-2±4.E-3 

a 20 

5.A07E-2+3.E-3 


t C2 
00 

2.405E+l±2.E+0 

a 40 

- \i 

3.566E-2i4.E-3 


,\2 

c fo 

-3.104E-2i4.E-3 


1.198E-3i4.E-5 


t 10 

-3.195E+0i3.E+0 

*£ 32 
r oo 

-5.374E-6il.E-6 

f? 

a ll 

-6.652E-2i4.E-3 

< 

3.197E-3i3.E-5 

* 

*31 

-9.5lSiS-3i2.E-3 

*Vo 

,7 

5.620E-3i3.E-4 


a 51 

1 . 114E-5i2.E-3 

*u 

5.582E-3+2.E-4 


1 

C li 

-2 . 238E-2±7 . E-3 

7? 

o 

O pj 

2.500E-3i5.E-5 


1 

C 21 

lo767E-2i4.E-3 

k a 
*20 

2.375E-3il.E-4 


b ll 

-8.004E-2+4.E-3 

1 

° c oo 

2.874E-2+3.E-3 


b 31 

6.322E-3±2.-E-3 

t cl 
00 

-2.726E+I±5.E+0 


b 51 

1 . 111E-2+2.E-3 

1 

C 10 

-l.660E-li8.E- 3 



-1.029E-3i7.E-3 

1 

C 30 

-2. 134E-2±8.E-3 


d 2 l 

-1.307E-2±4.E-3 

rt 

O O 
»— » 

-1.406E+0±2.E+0 

0 


a 22 

6.055E-4il.E-3 




b 22 

6.21lE-3i2.E-3 


*MSIS model value 






Table 2. Spherical harmonic coefficients for the atomic oxygen 
and helium concentrations at 120 km 


Coefficients 

0 


He 

Average (m”3) 

7,875E+i6i5.E+14 


1.981E+13il.E+ll 

a 10 

-l.178E-2il.E-2 

0 

(? 

-1.641E- 1 il.E- 2 

*20 <=> 

-3.763E- 1 il.E- 2 


-1.825E- 1 il.E- 2 

a 40 

-2. 18 IE- 1 il.E- 2 


>-2.278E- 1 il.E- 2 

f oo 

9.588E- 4 i2.E- 4 


-6.927E- 4 i2.E- 4 

*£ a ^’ 
_ 00 

0 


0 

if 

2. 26 IE- 3i2.E- 4 

X 

0 


1.892E-3 +8.E-4 


1.373E- 4 i3.E- 4 

*J a 
E ll 

2.160E- 3 il.E- 3 


6.400E- 3 i5.E- 4 

k a 

K oo 

-l.773E-3il.E-3 


-2.382E- 3 i4.E-4 

k a 
*20 

-8.695E- 3 +3 .E- 3 


-1.470E- 2 +8.E- 4 

•«S 

1.178E- 1 il.E- 2 


-3.647E- 1 il.E- 2 

o 

rt 

O 

O o 

b-* 

2.157E+1 i4.E+0 


-1.565E+ 2 il.E+0 

1 

C 10 

3.430E- 1 i2.E- 2 

l 


1.730E+ 0 ±2.E- 2 

1 

c 30 

1.272E- 2 ±2.E- 2 


1.933E- .1 12. E- 2 

cl * 
C 10 

5.949E+0 i3.E+0 


-5.945E+0 +5.E+0 

2 

c 00 

1.251E- 1 ±7 . E- 3 

■ (1 


4.093E- 1 17. E- 3 

2 

C 20 

-1.268E- 1 il.E- 2 

' « 

-2.073E- 1 il.E- 2 

t c2 
00 

1.030E+ 2 i2.E+ 0 


1.069E+ 2 17. E- 1 


O 


Table 2, (Continued) 
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$ 

Coefficients 0 He 


2 

10 

-l.574E-lll.E- 2 

-8.549E-2li.E-2 

CM 

a o 

i-l 

-8.260E+ 1 12.E+0 

-6.594E+U5.E+0 

l ll 

-1.273E- 2 tl.E- 2 

-7.072E-2i2.E-2 

l 31 

-1.826E- 2 16.E- 3 

-5.019E- 2i8.E- 3 

l 51 

-l.163E-2l4.E-3 

-5.425E-2l6.E-3 

1 

11 

1.813E- 2 12.E- 2 

2.380E-2 + 3.E- 2 

1 

21 

2.223E- 2 ±1 ,E- 2 

O \} ■ Q 

4.914E- 212.E-2 

'll 

3.600E- 2 ll.E- 2 

4.618E- U2.E-2 

31 

-6.465E-2 16%- 3 

-3.722E- 2 ±9.E-3 

51 

-3.277E-2 14.E-3 

-4.316E-2i6.E-3 

'll 

5.734E- 212.E- 2 

1.319E-113.E-2 

l 21 

2.339E- 2 li.E- 2 

2.220E- 112.E-2 

l 22 o 

8.638E- 3 ±4.E- 3 

2.404E-216.E-3 

22 

-3.489E-3 15.E- 3 

2.181E- 3 ±7.E- 3 


*MSIS Model Value 
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for consistency* These coefficients were determined in the S'jme way as 

\\ 

those for the exospheric temperatures* First, a quiet model was formed 
using 4*tu taken when the Ap index was less than or equal to ten* Then, 
the magnetic activity Coefficients were determined using all data with Ap 
< 130. No measurements obtained below 190 km were used. The reason for 
this is that atomic oxygen and helium very likely show departures from 
diffusive equilibrium at low altitudes* Since the model represents the 

ff 

hydrostatic distributions, a correction term for a non-zero diffusion 
velocity, such as that used by Hedin et a!. (1977b), would be needed in 
order to incorporate low altitude measurements into this model* No 
attempt to use a non-zero diffusion Velocity was made here* The 120 km 
lower boundary was chosen because it is a convenient level for modeling 

ft 

the upper thermosphere through the Walker (19(55) profiles. Actual 120 km 
values may be above or below the model estimates* 


As in the case of molecular nitrogen, the atomic oxygen and \hellum 
data were limited to low solar activity conditions so that the effect of 
EUV variations could not be found simply. In order to allow the model to 
be used under higher solar activities, the coefficients of the MSIS model 
describing the density variations with solar flux changes were used here. 
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CHAPTER V 

* 

DISCUSSION OF THE RESULTS OF THE MODEL 

o <= 

cf 5*1 Accuracy of the Magnetic Model 

The overall performance of the magnetic coordinate model was examined 

by computing the residuals (data to model ratios) of the atomic oxygen atv6 

\\ ' ’■ 

helium measurements as functions of day number, magnetic latitude, 
magnetic local time, nnd magnetic activity (Ap) in order to search for 
systematic removal of density variations with these parameters by the 
harmonic expansion* No attempt to estimate possible correlations between 

parameters was made, A subset (one of five points) of the quiet time data 

\ \ 

(Ap<, 10) comprised the density ratio plots, except for the magnetic 
activity residuals which were computed using data with Ap values as large 
as thirty*.. Only measurements obtained between 200 km and 400 km were 
included since the satellite coverage was best over this height range* 

The atomic oxygen residuals for 1974 and the first half of 1975 are 
shown in Figure 15a. Most of the data lie within a factor of 1.3 of the 
model, although several points are more than a factor of two below unity* 
These could be due to the slow return of the densities to quiet time 
values after a magnetic storm while the Ap index returns much more rapidly 
as is discussed in Chapter VI. There does not appear to be a systematic 
change in the residuals with day number indicating that the annual and 

<0 ■■■ 

semiannual effects have been largely included in the model* 
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Magnetic latitude behavior of the residuals is plotted in Figure 15b. 

i * 

The average scatter is fairly uniform over the entire globe, although the 
northern hemisphere contains most of the data w&Lch differ substantially 
from the model predictions. 

i Density ratios for various magnetic/ local times appear In Figure 15c. 

rp “ ' ° 

There does not seem to be any systematic change in the average scatter 

O 

with MLT but thn overall deviations are smallest between 1200 MLT and 1800 
MLT and largest between 2200 MLT and 0400 MLT. The latter region is one 
in which steep density gradients exist resulting from auroral heating. 
Because of the sensitivity of, the 0 densities in this region to the 
magnetic activity level rapid fluctuations of the densities would be 
expected in this area. „ aV 


The scatter of the residuals for Ap values up to thirty (Figure 15d) 

° " : . 

does not increase greatly with magnetic activity which is somewhat 
V’ 1 surprising since this model was developed for primarily quiet magnetic 

conditions. There are more data which fall below a factor of 1.3 of the 
model for Ap values between ten and twenty but most of these lie les3 than 

o •' ° 

a factor of two from unity. Above twenty the number of points is much 
smaller so that fej/er points having large deviations are expected. 


For the majority of the data, the model is expected. to reproduce the 


data to within a factor of 1.3. Some points fall well outside this range 


and these may be related to the delay of the composition in response to 

rtf If 

. ,- J - O ■ . * ■ ■ 

magnetic heating 'which is noCj represented by the Ap index. 


6 
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The helium density retlos over the period spanned by the AE-C 
database (Figure 16a) show somewhat more average scatter than those 
for atomic oxygen but the number of points differing greatly from the 
mode"^ Is ov»ch smaller. The majority of the data fall within a factor 
of 1.5 of the model. The slight rise in the average estimate with day 
number may be due to the adopted solar activity coefficients or may be 
spurious , 

Average scatter of the density ratio with magnetic latitude (Figure 

16b) is greatest in the polar region. This is expected since the 

(( 

densities are the most variable in the auroral zones. 

'SS. 

The 0000 MLT to 0800 MLT region shows the largest deviations of 
the data from the model although several points near 1800 MLT are also 
much lower than estimated (Figure 16c) . These points lie near the 
high latitude minimum observed in the long-term density study, an area 
of very large density changes. The composition is very sensitive to 
magnetic activity here and the same problems occur in modeling helium as 
occured with atomic oxygen. 

£> 

The helium density ratios for higher Ap values are not appreciably 

further from the model in general than those for lower Ap values although 

Cl 5 

several points are well outside a factor of two from unity for Ap values 
above ten (Figure L6d) . These points are few in number, though, and may 
be caused by local effects. 


.1 
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The helium modsl represents the major portion of the da til to within a 


factor of 1.5 up to Ap valuta of thirty* 'fh* largest dlscVepsneies appear 


to be in the polar regions, mostly betwiaft 0000 HtT and 08C0 MLT, an area 
of strong magnetic heating associated with the westward ejocfcrojet. The 


power input into the region is relatively large as shown by Figure 17 
(Nisbet at al., 1978) and substantial chants in the composition during 
magnetic storms occur too rapidly to be adequately represented by the 
nodal in this region* ° J 


5*2 Comparison of the Magnetic Modeland the MSIS Geographic Modal 
with bonnltude Terms to the Data i 


The major differences between an empirical model organized in 


magnetic coordinates and ope developed in a,; geographic system might be 




expected to occur i,*"- / tt|e polar region since the greatest density 


variations appear 


magnetic heating 


in til 


ese areas as a result of restricted but intense 


n thi\ auroral zones. The two coordinate systems also 


\ 


differ the most dk high latitudes. If the density variations are not 
adequately descried by magnetic latitude and magnetic local time 
components, but manifest significant universal time dependencies, the 
data should exhibit^ periodic changes with universal time. No such 
periodic structure is apparent in either the atomic oxygen or the 


helium residuals above 60° N or below 60° S geographic latitude as 


shown in Figures 18 and 19. The MSIS model with longitude terms (Hedln 


et al., 1979) shows no significant periodicities either, of course, but 

P o 


the scatter of helium is somewhat greater. This may be the result of 


C/ 


using a larger database to construct the MSIS model. 




, . .1 Al -r - 2 




* UNIVERSAL TIME (HOURS) 

Figure 18 Ratios of magnetic model predictions (top) 
and estimates of the MS IS geographic model 
with longitude terms (bottom) to atomic 
oxygen data between 200 and 400 km for Ap 
values less than or equal to ten above |LAT 




UNIVERSAL TIME (HOURS) 

Figure 19 Ratios of magnetic model predictions (top) 
and estimate of the MSIS geographic model 
with longitude terms (bottom) to helium data 
between 200 and 400 km for Ap values Igss 
than or equal to ten above |iLAT| of 60 . 
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For the AE-C data, the magnetic model seems to reproduce the high 


latitude composition variations at least as well as the HSIS model with 


longitude terms while using fewer parameters* Zt Is also much more 


convenient to use In particle and Joule heating studies since density 


estimates can easily be made In the coordinate system to whic^i these 


energy deposition processes are related* This is particularly important 


for comparison to theoretical models, where the number of coordinates Is 


limited i 


The low latitude behavior of the magnetic model is also quite 


important and must be considered since the major energy input in the 


equatorial region, solar heating, is geographically conS^lled. This 


energy source has a much more uniform global structure than the auroral 


heat input and the main variations can be represented using lower order 


harmonica. For these terms the magnetic and the geographic coordinates 


are about the same towards the equator so that either system should be 


satisfactory for modeling the composition. The density ratios of data 


taken within 30 of the geographic equator for the magnetic model show no 


significant periodicities in UT (Figures 20 and 21) and, in fact, are 


generally smaller than the residuals associated with the HSIS model, 


possibly for the same reason that the high latitude residuals differ. 


From the scatter in these residuals, it can be seen that the magnetic 


model represents the density variations as well as the MSIS geographic 


model with longitude terms at low latitudes. 


— -- . * . . 




UNIVERSAL TIME (HOURS) 

Figure 20 Ratios of magnetic model predictions (top) 
and estimates of the MS IS geographic model 
with longitude terms (bottom) to atomic 
oxygen data between 200 and 400 km for Ap 
.values less than or equal to ten below |LAT 
of 30°. 
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5. 3 Predictions of tha Magnetic Modal, the MSfS Geographic Modal 
without Longitude Terra*, and the MSIS Geographic Modal with 
Longitude Terms 


A comparison between the magnetic model estimates and those of the 

i) 

two geographic models by contour plotting the 120 km distributions in 
magnetic latltuda~uiagnetlc local tXuo; coordinates illustrates clearly the 

V h— ■ 

major differences between the moduli*-/ 


I 


Figures 22 , 23, and 24 show the helium densities for all three models 
for the same solar And magnetic activity conditions in the magnetic 
coordinate system* the predictions of tha MSIS geographic models for 


Various longitudes and universal times were, transferred to the magnetic 
coordinate system and put into bins of one hour MLT and 5° THAT and then 

. ’ - u O 

averaged in order to make these plots* The high latitude structure is 
quite similar for the magnetic model and the MSIS geographic model with 
longitude terms, but it is almost completely smoothed out in Jtha MSIS 
geographic model without longitude terms. The importance of magnetic 
position in modeling the neutral concentrations in the polar regions is 
seriously underestimated in the geographic model which has no longitude or 
universal time dependencies* 

o 

Although the MSIS model with longitude terms has the capability of 
reproducing both universal time and longitude variations simultaneously, 

i) ’ 

the magnetic model reproduces the main features of the data with fewer 
terms* It also does not smooth out the high latitude density variations 
as does the MSIS model without longitude terms, which has about the same 

cp .. 

number hi coefficients as the magnetic model* The advantage of having a 
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MAGNETIC LOCAL TIME 


Figure 23 Contour representation of log^o n(He) at 
120 km In the southern hemisphere for the 
MS1S geographic model with longitude terms 
at spring equinox for Fjq y • 100 and Ap • 








MAGNETIC LOCAL TIME 


Figure 24 Contour representation of log l0 n(He) at 
120 km In the southern hemisphere for the 
MS IS geographic model without longitude 
terms at spring equinox for Fiq 7 • 100 
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CHAPTER VI 


EFFECTS OF MAGNETIC ACTIVITY OH THE NEUTRAL COMPOSITION 

o 

6* l Resnonae of Atonic Oxyaen to Enhanced Man no tic 
Activity In the Rtaion of the Neat ward Auroral 
Elcctrolat 

Joule heating In the regions of the auroral electrojets increases 
markedly during magnetic disturbances* Fluctuations in nag nc tic activity 
throughout these disturbances can be monitored by the auroral electro jet 
indices (AE, AL, and AU), which attempt to measure the B field pertur- 
bations produced by the electrojet currents and by the 3-hour ap or 
daily Ap '’indices , which reflect global changes in magnetic activity, so 
that changes in the neutral atmosphere with magnetic activity variations 
may be modeled. Direct use of these Indices in modeling has, however, 
several problems. There is a time delay between the Joule heating and the 
winds and density perturbations and a further time delay before these 
affect densities at lower latitudes, and this results in time delays 
and storage phenomena (Wydra, 1975; Hedin et al. , 1977c). Therefore the 
temperature and composition of the neutral thermosphere cannot be repre- 
sented on a one-to-one basis with changes in the magnetic activity 
indices. Tills is particularly important in the polar regions where the 
energy input from magnetic storms is large and the density variations are 
substantial. Large errors occur in these areas in empirical models 
because data taken when the magnetic activity level is low immediately 
after a storm are used to develop quiet time models without the time 
history of the atmosphere being considered . 
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Thia problem la illustrated in Figure 25 whore the atonic oxygen 

o 

denalty at 120 km and the 3 -hour ap index have been plotted through the 
magnetically disturbed period 16 March, 1974 to 21 March , 1974, The 0 
densities were determined in the same manner as in the long-term analysis 
and were restricted to 55 to 60° ttLAT and 0400 to 0700 HIT In order to 
reduce density variations with magnetic position* Although the 3-hour ap 
index indicates the average trendc in the density variations, the 
relations with the 0 density during the recovery from the disturbance of 
16 March la very poor* The difference between the 0 densities here and 

those several days later after the atmosphere has recovered from the 

\ 

magnetic disturbance is greater than a factor of two, but no indication of 

thia slow recovery is given by the 3-hour ap index which is nearly* 

% - " ' ■ -V 

constant over this period* 

The same type of behavior was found in the recovery phase of another 
magnetic storm* Figure 26 shows the 0 density at 120 km between 55 and 
60° MLAT and 0000 and 0300 MLT and the 3-hour ap index for the period 10 
April, 1974 to 12 April, 1974* The densities do not return to quiet time 

values for approximately 24 hours after the 3-hour ap shows low activity* 

IP 

Again, there is more than a factor of two difference In the 0 values over 
this period, but essentially no change in the 3-hour ap index. Because 
such wide density variations with almost no change in ap values are 
responsible for very large errors in the empirical models, it would be 
desirable to include the atmospheric storage phenomenon in a magnetic 
activity indicator which could be used in modeling so that better, 
estimates of the neutral densities throughout and during recovery from 



UNIVERSAL TIME (FROM 0 HOURS, 3,6 MARCH, 1974) 

Figure 25 Logj^g n(O) at 120 km and 3-hour ap index as a function of universal time from 
16 March 1974. . 








magnetic » tor me could be determined* In tho region of Che westward 
electrojet, the AL Index la most nearly related to the ionospheric current 
system responsible for Joule heating and consequent density changes in the 
neutral atmosphere* Therefore, this index will be used exclusively in 
developing a magnetic activity indicator which will be related to atomic 
oxygen density changes near the westward electrojet for the two periods 
examined previously* 

ii 

i, 

6*2 Halation of a Magnetic Activity Indicator Based on 

the Auroral Electrolet Index AL to the Atomic Oxygen . 

Penalty Purina; Magnetic Dlsturbancaa 

An indicator of magnetic activity which is closely related to the 
neutral density variations which occur during magnetic atoms must include 
both tine delay and energy storage effects* At high latitudes, the delay 
between a magnetic disturbance and the resulting changes in the neutral 
atmosphere appears to be on the order of one hour (Prolss and Pricke, 
1976)* To a first approximation of a magnetic storm lasting more than a 

)'/ . o 

/? 

few hours, this delay time can be neglected. 

Storage effects can be approximated by using the AL index in an 

( 6 . 1 ) 

where the time constant, r , can be adjusted to obtain closest agreement 
between this magnetic activity index and the neutral density variations* 


equation of the form 


!lL(t) t « i J ^ AL<t')«~ (C " t#) / T dt' 


If 


The atomic oxygen densities at 120 km for the period 16 March* 1974 
to 21 March* 1974 and this index for t« 1 hour are shown in Figure 27* The 
index responds aa quickly as the atomic oxygen densities to magnetic 
activity increases but returns to qulst time values much more rapidly than 
n(0) Increases after the disturbance ( suggesting a value of t which is 

' 1 G 

much too small* For t* 12 hours (Figure 28)* the agreement between the Ml* 

index and the 0 densities at storm onset remains good and the relation 

6 

during recovery is much closer* although the densities still Increase more 
slowly than the ML index decreases after the disturbance has passed* A 
value of 24 hours for x (Figure 29) causes some delay in the ML index to 
the onset of the storm, but gives an even, closer relation between the 
index and the atomic oxygen densities during recovery. This relation is 
significantly better than that obtained using the 3-hour ap index to model 
the density variations . 

The ML index for t"12 hours also agrees more nearly with the 0 
densities than did the 3-hour ap index during the recovery phase of the 
storm of 10 April, 1974 to 12 April, 1974 (Figure 30) although the ML 
index returns to quiet time values more rapidly than they© densities 
increase* When T is lengthened to 24 hours, however, the atomic oxygen 
density values increase more quickly than-' the ML index decreases during 
recovery (Figure 31). An upper limit of less than 24 hours for Tin the 
ML index is thus Indicated for this storm* This value la consistent with 
chat obtained for the recovery phase of the first storm, though it appears 
to be larger than required to represent the onset. A value of t not 
exceeding 12 hours appears to characterize this phase more closely* 
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CHAPTER VII 
CONCLUSIONS 


7* 1 Lontt- term Averages of Atomic Oxygen and Helium 

Averages of the atomic oxygen densities reduced to 120 km and the 
helium densities normalized to 300 km for each season have been examined 
in a magnetic coordinate System for various magnetle activity levels under 
both low and high solar activities* A persistent depletion region for 
both constituents has been Identified at high magnetic latitudes (>60°) in 
the postmldnight magnetic local time' sector for all seasons and all 
magnetic activity levels examined, except for helium In winter* It is 
believed that the large seasonal variation In helium, at least a 25: 1 
ratio of winter to summer densities at 300 km, obscures this prominent 
feature* This region of low densities seems to be a permanent feature of 
the polar ^thermosphere since It is present under both high and low solar 


activity''! conditions* 

The location of the depletion region near the westward electrojet 
supports the theory that Joule heating due to electric currents flowing in 
the ionosphere causes upwelling of the neutral gas and a density decrease 
in the region of heating for constituents having less than the mean 
molecular mass (Mayr and Volland, 1972c; Reber and Hays, 1973)* 




120 


Largo density gradients in magnetic latitude and magnetic local time 
near the area of low densities suggest that the localised energy sources 
are comparable to or even greater than solar UV and EUV Inputs in these 
regions. 

The relation of the density minima to magnetic coordinates indicates 
that this coordinate system is well suited for modeling the neutral polar 
thermosphere* 

7* 2 Global Empirical Density Model 

I 

An empirical model of the atomic oxygen and helium distributions 
above 120 km has been developed in magnetic coordinates and compared to 

o'. 

similar models based on the geographic coordinate system* 

The magnetic model and the MSXS geographic model with longitude terms 
make similar predictions over the globe and show nearly equivalent 

structure towards the magnetic poles* The MSXS geographic model without 

!/ ' 

longitude terms does not show the same structure as these models at high 
magnetic latitudes, but rather gives highly averaged density estimates* 

Although the MSXS model with longitude terms has the capability of 
reproducing both longitude and universal time variations of the densities 
simultaneously, the magnetic model predicts very similar density values 
while using fewer parameters* The magnetic model is also based on a 
coordinate system in which many Important high latitude energy Inputs are 
organized which makes it a very convenient tool for use in a wide variety 
of studies of the polar thermosphere* J 

... ••••■- 

\ ■ 




A 1 










121 


7.3 Auroral Activity and the Atonic Oxygen Distribution \ 

r \ 

A magnetic activity Index baaed on the auroral electrojet lndek Al 

was developed for comparison with the atomic oxygen deta taken during two 

' \ 

magnetic storms. The auroral electrojet Index la a measure of the high 

\ 

latitude Ionospheric current and hence of the Joule heating In this 
region. The magnetic activity index was of the form* n 


MHt) - if* AL(t')e“ (t ‘ t,)/T dt' 

T T J -OQ f 


(7.1) 


Values of the time constant, t * were found to range from 12 to 24 
hours for storm recovery but leas than 12 hours for storm onset. 


7.4 Suggestions for Future Research 

One of the most important areas where further work In global 
empirical modeling Is needed lies in the response of the thermosphere to 

magnetic disturbances. Current magnetic activity indices such as ap and 

; ) 

Kp are not well correlated with the energy input to the auroral regions 
and to the resultant Joule heating. Major improvements in magnetic 
activity corrections must be based on indices more directly related to the 
actual energy deposition processed. The best current indices are the 
auroral electro jet indices AU, AL, and AE. Work should be done on the 
theoretical and experimental relations of these indices to the spatial and 

temporal variations of the energy inputs and on the atmospheric response 

0 ' 0 

(Niabet, 1980). 


0 


Measurements of lower altitude thermospheric temperatures oyer the 
entire globe which are direct would also be of tremendous value* Such 

information would permit variable lower boundary conditions baaed on 

o 

actual data to be employed in the models and, thereby/ expand enormously 
the altitude range over which model predictions would be realistic. 

Density and wind velocity data in the same region would be 
instrumental in determining the effects of turbulence on the composition 
and would serve as a valuable link between photochemical equilibrium 

o 

models at lower heights and diffusive equilibrium models at satellite 
altitudes* 


!) 
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